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(54) Methods of characterising ventricular operation and applications thereof 



(57) New methods of characterising ventricular op' 
erations by measuring pro(3agation characterisics of the 
repolarisation wavef ront (the T wave) are disclosed the 
methods use new descriptions of T wave Morphology 
Dispersion (TMD), Total Coslh R_ to_T (TCRT) and T 
wave energy res idium to quantify the wavef ront charac- 
teristics, these descriptors measure the spatial variabil- 



ity of the T wave Morphology, the vector deviations be- 
twe.en the depplarisation and repolarisation wavef ronts 
and the energy of the non-dipolar cdmpbnerits of the 
ECGvector respectively TCRTaIsp provides a respon- 
sive descriptor for measuring autpnorriictone As sucii, 
has applications for;, improved pacing and autonomic 
nenrous system monitors. 
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Description 

[0001] The present inventions relate to methods of characterising ventncular operation In particular, but not exclu- 
sively, they relate to a system for quantrfymg abnormalities of an electrocardiogram and to a method and an apparatus 
5 for measuring such abnormalities. The present inventions also extend to an operaiing system for a computer, to a 
computer program and to media having stored thereon a computer program for putting the inventions into effect. Other 
applications include, use of the algorithms in pacemakers and heart monitors. The inventlons share a common lirik of 
characterising differences in the wavefront of the repolartsatiori wave. 

[0002] Electrocardiographic patterns of the heart's movements have been well studied. An electrocardiogram (EGG) 
10 records the changes in electrical potential associated with the spread of depolarisation and repolarisation through the 
heart muscle in a normal healthy patient, depolarisation starts in an area of the right atrium called the sinoatrial node 
and spreads through the atrioventricular node and into the ventncular muscle via specialised conduction tissue, causing 
the two atria and the two ventricles to contract Dunng repolarisation, the atria and ventricles relax and refill with blood 
The depolarisation of the atria is responsible for the P wave of an EGG and depdiarisatipn of the ventricles results in 
IS the QRS complex Repolarisation of the atria coincides with the QRS complex so it Is not seen Repolarisation of the 
ventricles, however, is seen as the T wave 

[0003] EGG'S are typically recorded using a standard arrangement of 12 leads, 6 (the I, N, l|l, VR, VL, VF leads) 
looking at the heart In different directions in an approximately vertical plane of a body in an upright position and 6 (the 
VI, V2, V3, V4, V5 and V6 leads) looking at the heart in different directions in an approximately horizontal plane Using 
20 such an arrangement of leads, the spread of the waves of electrical potential associated With depolarisation and re- 
polarisation through the three dimensional space of the Ixady, can be retSorded. 

[0004] The spread of these waves through the heart is often described by vectors. For example, the average direction 
of spread of the depolarisation wave through the ventricles as seen from the front of the body Is.called the cardiac axis 
and Ihe direction of this axis has long been used to indicate different abnormalities of the heart 

2S [0006] To study abriomnalities associated with ventricular repolarisation, a number of data processing techniques 
have been proposed to measure, for exannpie, the QT interval, i s the ihtewal between the beginning of depolarisation 
and the end of repolarisation of the ventricles Interlead variability of the QT interval durations in standard 1.2 lead EGG 
recordings has also been studied However, whilst these measurements may provide some diagnostic assistance, 
concerns have been raised about the poor reproducibility of resuKs 

30 [OOOS] Studies have also tried to quantify the inhomogeneiiies in the ventricular repolarisation patterns by ev9lu9ling 
the complexity of the T wave morphology using eigenvalues associated with the principal components of EGG, nneas- 
ured over a period of 24 hours The direction of the EGG vector during T wave in the 3D physical (x,y,z) has also been 
shown to have gome predictive value 

[0007] However, there is still a need for further measurements which may provide a more accurate technique for 
3S identifying certain conditions, particularly those which affect repolarisatiw of the ventricles A problem with known 
methods, for example, is that they only quantify global variations in the T wave rather than spatial variations in indiwdual 
waves, that is the synchronicity of tie T wave, as observed from, different locations on the body, is not observed 
[0003] Thus, viewed frorrt a first broad aspect, a first inventbn described herein provides a method of quantifying 
abnonnalities of an electrocardiogram observing repolarisation patterns from different locations on a body, wherein the 
40 abnormalities are quantified by a measure of the synchronicity of the repolarisation patterns as. observed from those 
different locations on the body In other words, this Is a measure of the homogenily of the spread of repolarisation waves. 
[0009] Unlike depolarisation of the heart m jscie, the repolarisation of individual cells is not triggered by neighbouring 
cells but IS instead a time dependent process If repolarisation patterns, as observed in different locations on the body, 
lack synchronicity, then this can be Indicative of certain heart complications 
46 [001 0] By quantifying these abnormalities, it may be possible to use the data to assist with diagnosis or to identify 
patients at most risk or classify them Into different categories of risk. This may be of great importance in determining 
whether certain treatments should be offered to a patient, for example The data could also be used to trigger an alert 
in a monitoring device 

[0011] The homogenlty of the spread of repolarisation waves can be measured by quantifying the spatial variability 
so of the ventncular repolarisation patterns i e the spatial variability of the T wave 

[0012] Thus viewed from a second aspect, the first invention provkJes a method of characterising ventricular oper- 
ation, comprising the steps of: 

recording a signal monitoring the propagation of a repolarisation wave, 
ss determining a vector which is representative of the wavefront of the. repolarisation wave, and 

determining a measure of the spatial variation of the repolarisation wavefront. 

[0013] in one preferred embodiment, it providesd a method of quantifying abnormalities of ventncular repolarisation 
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by determining a measure of thG spatial T wave morphology variation; 

[6014] Preferably the spatial T wave morphology variation is quantified by measuring the T wave Morphology Dis- 
persion 

[001 S] Preferably this is achieved by detemiining vectors describing the contributions which the signals from each 

s lead (often referred to as the channels of an ECG) makes to the T wave. The angles between these vectors are then 
calculated and a mean value is determined: This mean \^lue of the angles provides a measure of tfie spatial T wave 
morphology variation. The smaller the value, the closer the T wave morphologies will be in the signals of the individual 
■ leadS: ■ " -■ : , 

[0016] Preferably the ECG signal ismorphologically filtered to improve the signal to noise ratio, in one preferred 

10 embodiment, this consists of the steps of decomposing the T wave using a technique such as Singular \Ja\ue Decom- 
position filtering by keeping only the two most significant signal components, and applying a DC compensation A 
preferred DC compensation is provided by. subtracting an average of the start and end signal components during the 
QRS complex and T wave The rnbrphoi.6gically filtered T wave is then preferably rescalsd to equalise energies in the 
different component dirBctbns. The corresponding reconstmctipn parameters are calculated to determine the vector 

15 contributiOTts of each of the ECG leads. The angles between each pair bf the vector contributions is tiien calculated 
and Ihe mean determined Most preferably the contribution in respect of lead vi is ignored because the T wave mor- 
phology In this lead is generally different than that of other channels, Irrespective of any clinical background, mainly 
due to the position of the VI electrode, and by ignoring this component, it has the effect of enhancing tfie predictive 
value of the T wave morphology dispersion descriptor ■ 

:S0 [0017] The main reason for initially decomposing the data matrix is to find an optimum representation of the ECG 
signals upon which the measurements can be performed In this way, the system does not use the standard XYZaxes 
of the body but finds an optimally constructed orthogonal system to representthe 12 lead ECG In a preferred embod- 
iment, therefore, the first Invention can be seen as providing a method for .iooking..at the yectpr representatioih pf each 
of the standard electrocardiographic leads in an optimum dimensional, vector space in which the ECG signaSs pan be 

25 represented and comparing the angles between the vectors Gfindividual standard leads. . . 

[0018] The spatial T wave. morphology variation may provide a useful descriptor w^en It is determined for the whole 
of the T vkve, the first half of the T wave, the second half of the T wave or any other portion or combination of portions 
of the T wave 

[0019] The present inventions are not limited to standard 12 lead electrocardiograms, although this is preferred, but 
30 extend to electrocardiograms produced from only throe or more leads In certain applicalions, it may be useful to use 
the electrodes of a pacemaker to record an electrocardiogram signal In such situations, the positions and numbferS of 
the electrodes would hot usually correspond with the arrangement of standard leads, for example Whilst It is preferred 
to view the waves in three dimensions, because, research upto-date suggests that iapproximately 99% Of 1 2-|ead ECG 
energy can be represented Ih a 3D space, the Inverrtions are applicable to situations where the .heart is viewed in any 
35 dimensions, greater than or equal to two, 

[OMO] In a conventional 12 lead ECG, only 8 (I, II, VI , V2, V3, V4, V5 and V6) Of the signals are independent The 
Other 4 signals (III, VR, VL, VF> are algebraically dependent on the other leads so, if desired, may be generated by 
data processing methods rather than measured as such As explained above, it is most preferred to use signals only 
from leads I, I!, V2, V3, V4, V5 and V6, and to ignore the signal trom VI in order to concentrate the abnormalities seen 
40 In the T wave' The position of the leads, although having an effect on the value produced by the descriptors. Is not 
critical to the inventions. Whilst the inventions may bo described with reference to the standard ECG leads, this is not 
intended to limit the inventions to just those positions These positions are preferred, however, since they have become 
"standard" measuring points throughout the world 

[0021] It has also been found that comparing the spread of dspolarisation through the ventncles with the spread of 
45 repolartsation can provide usefufinformatlori 

[0022] Thus, in accordance with a firist aspect of a second Invention described herein, there is provided a method of 
characterising ventricular qperation, corhprislng the steps of; 

recording a signal monitoring the propagation of depoiarlsation and repolarisatiori waves, 
so determining vectors which are representative of the direction of the wavefronts of the depplarisation and repolari- 

satfon waves, and 

determining a measure of the deviation between those vectors 

[0023] Thus, the present invention can be seen to provide a me»iod of quantifying abnormalities of an electrocardi- 
55 ogram observing the spread of depolarisation and repolarisaticn waves through the ventrteles, wrfierein the abnormal- 
ities are quantified by comparing a property of the depoiansation wave with a property of the repqiarisation wave where 
preferably the abnormalities are quantified by a measure of the vector deviation between the ventriculardepolarisation 
and the ventricular repolarfsation waves 
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[0024] Described in other terms, a method of the second invention may compare the direction of the depolarisation 
wave (1 e , the QRS part of the ECG) with the repolarisation wave (i e the T wave) This may be achieved by comparing 
the angles between principal vectors of the ventricular depolarisation and repolarisation waves, comparing. tie angles 
between a principal vector of the ventricular depolarisation wave and the ECG vectors during ventricular repolarisation, 

s connparing the angles between the ECG vectors during ventricular depolarisation and a principal vector of the ventricular 
repolarisation wave, or comparing the angled between the ECG vectors during ventricular depolarisation with those 
during ventricular repolarisation. The angles may be comparesd for the whole of a wave or just a portion of a wave or 
any combination of portions of the waves: For example, in one embodiment the angles between the depolarisation and 
repolarisation vectors are compared for por1;iqns of the waves which span the peak energy values, but it may be pre- 

10 f erred in some instances to look at and opmpare other pdrtbns of the waves which would correspond to depolarisation 
and repolarisation occurring in different regions of the heart muscle. 

[0025] Preferably the vectc^ deviations are determined using the optimally constructed representation of the ECG 
signals discussed above 

[0026] In healthy patients, the principal vectors would, generally speaking, only deviate by up to about 30" , In patients 

T5 having hypertropic cardiomyopathy (HCM), for example, vector deviations greater than 90° may be seen These de- 
viations can be distinguished over inversion of the T wave, for example, which would result in angles closer to 180° 
[0027] Again it is preferred to conduct certain data processing steps before the angles of the vectors describing the 
QRS complex and T wave are compared Firstly the data matrix describing the signal Is decomposed, again preferably 
by Singular Value Decompositim The decomposed signal conporients are ranked In order of their significance In 
terms of the energy of the ECG vector that they represent Thus the first signal component contains the most energy 
in a first direction The second signal component contains the next most energy in a second direction which is perpen- 
dicularto the first. The third signal component contains the next most energy in a third direction which is perpendicular 
to the first and second directions, Where eight of the independent ECG channels are recorded, the ECG vector can 
be decomposed into an eight dimensional orthogonal space When measuring the vector deviations.of depolarisation 

2S and repolarisation vectors, a good approximation may be made by only measuring the first two or three of the depom- 
posed signal components since these can account for more than 99% of the total energy of the 12-lead ECG signal 
The QRS complex aiid T wave are localised by making use of the variation of the instantaneous ECG energy. The 
method does not depend on accurate localisation of the QBS complex and the T wave This method of detection 
provides an example of many possible ways 

30 [0028] Both the visctor representations of the QRS complex and T wave follow an approximate loop, in the constructed 
space Vectors can he determined which represent the maximum energy of the T wave and QRS complex, and the 
angles between therh compared More preferably a vector describing the maximum energy of the T wave is compared 
to the vectors describing the QRS complex for a set of points at the peak of ventricular r^olarisation (which, for 
example, can be determined by tracing the instantaneous ECG energy), in the most prefen-ed embodiment, the vector 

3S deviation between the depolarisation and repotansation waves is measured as the average of the cosines of the angles 
between the vector describing the maximum energy of the T wave and the vector describing tlie wavef ront of the QRS 
complex, tlis angles being determined in corrstriicted space The measurement of this vector deviation In terms of the 
cosine of the angle is referred to herein as TCRT - total cqsine R_'to _T 

[0)29] The second invention introduces the idea of considering depolarisation and repolarisation processes of the 

40 heart muscle sim ultaneously and comparing them While it may be described as a comparison of ECG vectors observed 
during these processes, in an appropriately constructed vector space, this is not intended, to limit the invention to just 
ECG vector comparison and to processing of standard ECG leads 

[0030] Coriiparin.g the propagation directions of depolarisation and repolansation has revealed some interesting de- 
tection properties In particular TCRT has proved to be more sensitive to autonomic changes of ventricular repolari- 

4s sation than other known descnptors such as ventricular gradient and QT dispersion From investigations, it has been 
found that TCRT responds quickly to changes in Bie position and activity of the patlerit with distinct ranges or levels 
of descriptor values being obtainable for different autonomic tones This descriptor may be used to check for abnor- 
malrties in ventrcular depolarisation and repolarisation under different autonomic conditions, thereby providing a fuller 
picture to assist with diagnosing defects TCRT has been found to be useful as a predictor for mortality in patients 

50 which have suffered acute myocardial infarction and as a predictor for arrythmias. The descriptor could be used in a 
monitor carried by the patient or in equipment in an intensive care unit to warn the patient or medic by means of an 
alarm when the TCRT is pushed to a danger level as a result of changes in the autonomic tone, for example, caused 
through exercise or trauma. TCRT could also be used to check that the patient has a properly functioning autonomic 
system, for example, prior to the administration of anaesthetic agents before surgery. 

55 [0031] TCRT is in effect able to provide a measure of the autonomic tone of a patient It could be used to control 
pacing of a pacemaker making it more responsive to the patient's needs by responding to changes in the autonomic 
tone TCRT is increased by physical or emotional stress during fixed rate pacing and is decreased by an increase in 
pacing rate TCRT could be implanted in a closed loop rate-adaptive feedback system At times of physical or emoticml 
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stress TCRT would increase, triggering an increased pacing rate to decrease TCRT to resting level. 

[0032] TCRT could be used to monitor the effect of certain drugs and the way in which they effect the autcffiomic 
system When testing drugs which prolong the QT interval, TCRT could be used to monitor the patient and raise an 
alarm arrythmias are predicted or detected. It could be used to monitor changes in electrolyte of the body and other 
5 conduction phenonoma. If could even be used to control a drug delivery mechanism, administering certain drugs as 
heart f unction abnormalities are detected or in response to changes in the autonomic to^ 

[0033] Conditions such as ischemia, as well as most illnesses, will have' an effect on the autonomic tone of the 
: patient. TCRT could be used to assist in the prediction of ischemlaor in the monitoring of the progression of a disease, 
for example, in heart failure patients, by providing an indication of the autonomic tone as well as changes in direction 
10 of the re polarisation wavefront. TCRT may be able to obsen/e autonomic changes caused by the onset of ischemia 
before ST segment changes are observable on an ECG or pain is felt by the patient It may also be useful tn the 
monitoring of patients sufferingf rom epilepsy, providing an early warning of heart function abnormalities The autonomic 
tone could be observed to detect hypoxy conditions in a patient prone to fitting. 

[0034] . It should be noted that the possibilities mentioned above with reference to TCRT are not intended to be inferred 
IS as limiting the present invention to the preferred situation where TCRT is the average of the cosines of the angles 
between the vector describing the maximum energy of the T wave and the vector describing the wavefront of the QRS 
complex For example, TCRT may compare the angles between the sets of vectors describing the depolarisation and 
.t/repolarisation wavefronts with respect to time or may compare the vectors describing the maximum energy of depo- 
larisation and repolarisation to each other While the use of cosine, provides an effective way of separating the angles 
so . between the vectors associated with abnormalilies from those observed in normal patients, other operators may be 
used to separate the data 

[0035] As mentioned above, of the standard 1 2 lead ECG signals, .8 a,re independent Thus, It is possible tQ describe 
..the T wave as an 8-by-n matrix M, with each, row corresponding to a standard EC|3 channel (I, i\, VI, y2, V3. V4, V5, 
V6) and n being the number of samples Performing Singular \felue Decomposition on ftie matrix M generates a signal 

2S vector representing the progressof a T wave in S dimensions, where each dimension cari be regarded as a component 
of the signal, associated with a fraction of the total energy of 1 2-lead ECG For most purposes, as mentioned above, 
only the first two or three components are nonmaliy used, since these may account for over 99% of the total T wave 
energy In 12-lead ECG signals However it has been found that comparing the energy oJ the most Significant compo- 
nents with the energy of the other components provides a further useful descriptor that can be used during analysis 

SO of the ECG 

[003G] Thus according to a third invention described herein, there is provided a method of quantifying abnormalities 
of an electrocardiogram having a plurality of independent signals, in which the signals are decomposed to obtain a 
signal vector having two or more signal vector components, wherein the energy of the components is Compared. 
[0037] ' Preferably, the components are arranged substantially in order of decreasing signaf energy and the energy 

35 of the most significant components representing the majority at the signal energy is compared to the energy of the 
other components, Preferably :lhe electrocardiogram records 8 independent signals and the signal vector has 8 com- 
pcwients arid the first 3 components of the signal vector, thus constructed, represents the majority of the signal energy 
[0038] The third invention introduces the idea to transform the ECG signals into an optimally constructed space which 
represents the ECG energy in a minimum dimensional space of orthogonal components and to assess the residual 

40 energythat is left outside the three dimensional space of the first three signal components. In a possible implementation, 
as described above, this corresponds to assessing the relative values of the 3 highest singular values of the matrix Wt 
to the other singular values. The first three signal components represent the dipolar components of the ECG vector 
and the remaining signal components correspond to the non-dipolar components.. By measuring the power in the 
orthogonal components outside tfie three dimensional space, it is possible to measure the power of the non-dipolar 

45 components which provides a nrieasurs of the local repolarisaticsi abnormalities. However with a different set of Input 
signals and/or different parts o' the ECG signal, the separation/number of significant components may vary 
[003S] The inventions described above also extend to an apparatus which is programmed with an algorithm to proc- 
ess data from an ECG in accordance with any cif the described inventive methods The apparatus could be a computer, 
for example, programmed in a particular way or could be a plug-in box for an ECG apparatus or an ECG apparatus 

so provided with means to calculate these descriptors and oisplay the result to an operator Furthermore, the inventions 
extend to an operating system ora computer program having an algorithm to process data from an ECG in accordance 
with the described methods and to media having such a computer program or operating system stored thereon Thus, 
the inventions extend to a computer program product which is directly loadable in.to the internal memory of a digital 
computer, comprlstng software code portions for perfomilng the steps of the afore-dascribed methods when the product 

S5 IS run on a computer 

[0040] The present inventions wi II now be described by way of example only with reference to.a preferred embodiment 
and the accompanying drawings, in which. 
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Figure 1 shows a flow chart of praferred algorithms; 

Figures 2a and 2b shows 8 input and 8 decomposed ECG signals (which can be obtained by Singular Value 
Dacomposltion), respectively, 

Figure 3 illustrates approximate QRS and T wave detection using ECG energy. Ego, which is calculated Irom S,. 
5 and S3, the most significant 3 decomposed ECG signals; 

Figure.4prOTldes an. example of the T loop: Path of the tip of Sgo (the rndst significant point of the ECG vector in 

the deccsmposition space) on the u^Uj plane; 
V Figure 5 Illustrates an example of the cell weights D's, assigned to 100 equal cells in u^Ug plane (excluding the 

zero weights), in increasing order, showing the 30"> cell above the threshold, the 30* cell (the one shown in Figure 
10 4) is selected as the T wave end, 

Figure 6 shows the 3 most significant decomposed, time-orthogonal, channels of a T wave, 

Figures 7a artd 7b illustrate a T vrave loop and reconstruction vectors of each standard ECG lead for a normal and 

HO M patient respectively; 

Figures Sa and 8b illustrate QRS and T vyave loops for a normal and H.CM patient respectively, 
75 Figure 9 iOustrartes the reprodupibihties of all the descriptors, as measured by the ratio of individual variation to 

total variation in 1 0 supine recordings, from. 76 normal and 63 HCM subjects (Grey; Normal subjects. Black. HCMs), 

Figures 10a and 10b show the variations in descriptor value for ventricular gradient angle, ventricular gradient 

magnitude, TCRT, RR interval, maximum QT and CUT dispersion for the postural ch^ges of resting supine position, 

followed by sitting, unsupported standing, supine and standing position, 
20 Figures 11a and lib show the variations in descriptor value for ventricular gradient angle, ventricular gradient 

magnitude, TCRT, RR interval, maximum QT and QT dispersion during Valsalva manoeuvre, each manoeuvre 

being performed 3 times in supine position and 3 times in unsupported standing position, 

Figure 12 ^ows an implantable device positioned within the human body, 

Figure 1 3 shows a perspective view, of a preferred pacernaker, 
gs Figure 14 shows a circuit diagram suitable for operating the paceriiaker of Figure 13, 

Figure 15 shows a perspective view of a. preferred pacemaker-cardibverter-defibrillator; 

Figure 16 shows a circuit diagram suitaible for operating.ihe pacernakerof Figure 16, 

Figures ITaand 17b show scatter diagrams of QT dispersion values, where QTd Method 1 = range of measurable 
QT intervals, QTd Method 2 = standard deviation of measurable QT intervals, and QTd Method 3 = intenquartile 
30 difference of maasijrable QT intervals, 

Figures ISa and 1 8b show QTd Method 1 and QTd Method 2 results displayed with respect to the clinical group 
of the subjects, 

Figure t9 shows the relative T wave residuurn ptotted with respect to the clinical group of the subject, and 
Figures 20 and 20b show the lack of relationship between QT dispersion and relative T wave residuum data 

35 

[0041] Three new approaches forths analysis of ventricular repolarisation in 12 lead electrocardiograms (EGG) will 
now be discussed in relation to. a first case study: The spatial and the temporal variations of T wave morphology and 
the WavefronJ direction difference between the ventricular depolarisation and repolarisation waves A minimurr dimen- 
sional space, constructed by the Singular V&lue Decomposition of EGG signals, Is used The spatial variation charac- 

40 terises the moiphology differences between standard ieada The temporal variation measures the change of interlead 
relations during ventricular repolarisation The depolarisation and repolarisation patterns are compared using the prin- 
cipal Interlead relations (wavef ront directions) that charactenses them All of ttie descriptors are measured using the 
ECG vector in the constructed space and the singular vectors that define this space None of the descriptors requires 
time domain measurements (e g. the precise detection of the T wave offset), avoiding the inaccuracies associated with 

45 the conventional QT interval related parameters 

[0042] The new descriptors have been compared with the conventional measurements provided by a commercial 
system for an automatic evaluation of QT interval and QT dispersion In digitally recorded 12 lead ECGs (QT Guard, 
Marquette Medical Systems) The basic comparison used a set of 1100 normal ECGs The short-term intrasubject 
reproducibility of the new descriptors was compared to that of the conventional measurements in a set of 760 ECGs 

so recorded in 76 nonmal subjects and a set of 630 ECGs recorded in 53 patients with hyperthropic cardiomyopathy (1.0 
serial recordings in each subject of both these sets) The discriminative power of the new and conventional parameters 
to distinguish normal and abnormal repolarisation patterns was compared using the same set 
[0043] The results showed that the new parameters do not correlate with the oonventionai ones (i e. assess different 
ECG qualities), are generally more reproducible than the conventional paramaters based on the QT interval raeasute- 

55 ment, and lead to a more significant separation between normal and abnormal ECGs both univariately and In multi- 
variate regiession models 

[0044] The present Inventions are based on the hi^thesis that spatial and tennpora! variations in T wave morphol- 
ogies and the relatian between the depolarisation and repolarisation patterns will offer new measures of repolarisation 



6 



BNSDOCID: <EP 1 038498A2_I_> 



EP 1 038 498 A2 



abnormalities, it was aimed to define a set of parameters that would (a) quantify such abnormaiities. (b) be highly 
reproducible, (c) have sensitivity and specificity greater than the conventional measurements such as QTd in separating 
normal and clinically relevant abnormal electrocardiograms, and (d) be Independent of problematic time-domain meas- 
urements such as the detection of the T wave offset 

[0045] EGG decomposition by several different methods fs known. One example is taught by Acar B. and Koymen 
H In an article entitled "SVD-b&sed on-line exercise EGG signatorthbgonalization", (1 999) IEEE Trans: Blorned, Eng., 
46, pp: 311-321 . The common approach is to transform the multiple lead EGG into another domain or to identify the 
dominant components of the recording. The analysis is subsequently carried out in the transformed dorriain which may 

have various advantages, e.g.. a high signal-to-noise ratio.. - 
[0046] In the lollowing description, the uppercase bold letters are used for matrices whereas the lowercase bold 

letters for vectors 

[0047] The analysis is based on the Singular \felue Decomposition (SVD) of the standard 12-lead EGG It defines a 
minimum dimensional subspace that captures the EGG energy SVD was defined by Golub G H., and Van Loan C F 
(199S) in 'Matrix Computations', 3"^ edition, {The Jdin Hopkins University Press, Baltimore and London), pp, 70-71, as. 
[0048] If M is an B-by-n matrix, with each row corresponding to a standard EGG channel (I, II, VI , V2, V3, V4, V5, 
V6)and nbeing the number of samples, then there exist orthogonal matrices 

U = [u,,. . ,u„] eS?'""" and V = [v, , ...,v,} eyj""" 



where > og ^ ^ Op i 0 

[0049] The columns of U are referred to as the left singular vectors, whereas the columns of V are referred to as ths 
30 right singular vectors a, are the singular values. Furthemnore, if >. ^ o.^> ct^+i =•■•= «Jp. = 0, then 



■ range(M) = span{u,,.. ,u,} 

where rang^M) is the minimum dimehsioral spacewhich captures the 1 2-lead EGG energy As shown m the following 
section, no significant errors are introduced by restricting the dimension, of this space to 3. The singular values are 
measures of how much EGG energy exists along the corresponding vector u 

[OOSO] The signal representations in the minimum (r-dlmensional) decomposition sub-space and the corresponding 
left singular vectors arc used to derive the new parameters (Figure 1 ) in the rest of this text, iU designates the ECG 
data matrix. Each column of M corresponds to a sampling instant and each row corresponds to a different ECG channel 
Because of the algebraic interdependency, eight of th e staridard 1 2 ECG channels sire used, namely 1, 1 1 , VI , V2, V3, 
V4, V5 and V6 The SVD of lU is perfomned as. 



'to 













iVlV = U'"MV 

r 

£ e and is diagonal , U e g?"" 
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wheregtsi" and is diagonal, and It has been shown previously that:99% of the ECG energy can be represented 

In a 3-drmensional minimum subspace (see the paper mentioned above by Acar and Koymsn, 1999). Hence the ef- 
fective rank of M is 3 This mimmumsubspace is spanned by the columns of OSS'-' Let Sbetha projection of M onto 0 .. 
S =rTM Figure 2 shows an example of the input and decomposed (projected) signals. The transformation of the 3 
5 dominant decomposed signals back into the original ECG domain is equivalent to morjshologlcal filtering of the ECG 
in its original domain 

[0051] An approximate QRS and T wave detection Is perfoimed on the 3 (s, , s^, Sg) decomposed signals that contain 
most of the energy. Let : 

10 : 

[0052] Figure 3 shows s^Q anti Ejpfor a single beat The R-wave end point is assumed to be the first point after the 
maximum of Eg^^ where Ego falls below 70% (arbitrary threshold) of its maximum value That instant is disclosed as 

20 tfj^ Although fflipis not the actual R-wave end point, it serves the purpose Similarly, the 70% point before the maximum 
Is marked as fpg. The QRS complex is assumed to start 48 msec before f'R£and end 48 msec after t'f^^ (arbitrary 
limits) These two points are marked as f'^^and f'^jf, respectively The T wave peak tjp is assumed to be the maximum 
point of E30 after t'^. The approximate T wave starting point tjg Is taken to be 1/3 of (tjp - t'p,^ after f'^^. Figure 3 
shows the positions of t'^^, t'fj^, fpi^, f'^f, i^and tjp. Since this study Is only concerned with single beats, there isr ho 
need to choose an end point for the search of the T wave peak. However, such a point can easily be selected based 
on the instantaneous heart rate, when necessary The. detection of the T wave end point, fye is based on the path of 
the tip of Figure 4 shows the path that the tip of ^^{y follows in. the plane spanned by and Ug, for > (7-5. The 
rectangular area is divided into 100 equal rectangular cells (arbitrary number) Each cell is assigned a weight equal to 
the number of inner data points (time instants) This is a measure of the time spent by the tip of in that coti 1 fit D, 

30 be the weight of 1he i'^i cell The cells with zero weight ars discarded and the other cells are ordorod in respect o! Dj. 
Figure 5 shows Devalues for a single beat Assuming that there are cells with nonzero weights, < Do <-< 
[0053] A threshold th - mean{D,) + X standard_deviation(D,) is used, where ^i =3 (arbitrary constant) If i th 
forL<i<. Kior Y< i< O, the earliest time instant at which the tip of .Sgo enters one of the cells Vto O (L< i< K), is 
set to be the approximate T wave end point, tjB 

3S [00S4] .Since «ie aim of the algorithm is to quantify the T wave shape between tjs and tjE, rather than to measure 
the f'Rginten/al, the approximate and arbitrary nature of the r^s instant definitions is fully acceptable (as shown further) 
[0055] If, using the algorithm described above, tj^, which should not occur, |i is increased in steps of 0 2 until fy^ > 
tjp. Similarly, if < ih. ]i is decreased In steps of 0 2 until D/f > th. Such cases are rare 

[0056] This T wave end point detection scheme is based on the concept that the interlead relations do not change 
40 in the absence of the ECG signal Each point on the Uo-plane corresponds to a specific interlead relation defined by 
the vectors and Hence, each cell in the u^Ug-plane represents a group of similar interlead relations. When the 
repolarisation pattern ends, the ECG signal remains confined to a small set of such relations 
£0057] The decomposed signal is subsequently normalised with the maximum energy set to 1. 

45 

[0058] QRS complex and the T wave extracted by fhjs algorithm result in decomposed data matrices S^-and ^qbs' 
ss respectively From both signals, a DC vector is subtracted as folk>ws. 
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= 0.25 X {S,„(/,,,) + s,„ {t^^ ) + s,„iir^} + S,^ )} 



5 



[0059] Frcsn now on, 



S^ GSi^'^ {K'=t^--t„) and c 9?'"^ (i = /j,^ -/^j^r 



TO 



will denote the decomposed, energy normalisec and DC'Compensated T wave and QRS complex 
[0060] The T wave is reconstructed irom Sj, which is equ ivalent to a moiphological filtering: 



15 



tThe reconstructed T wave, M^-, Is once again decomposed by S\©: 



20 




IS diagonal ti^es-' ^^^x'" andisdiagonah.rSr^^'" 

[00S1] ThG subscript T indicales that we are dealing with the T wave only, and the superscript 'z^' denotes a matrix 
which is reconstructed (morphologically filtered) Note thato, has two columns whereas c has three columns This is 
because the S"' decomposed signal in T wave decomposition has been excluded ( Figu re 6) 

[0062] The spatial variation descriptors are determined as follows, t), is an 8-by-2 matrix Its columns are the two 
most sigmficant leftsinguiar vectors o1m„ 



Each of Its rows Is the reconstruction vector of the corresponding standard ECG lead (note that the reconstruction from 
the 2D, most significant, subspace of the decompositiwi ispace means multiplying the decomposed data matrix by ) 
Let z denote the reconstructbn vectors, 



45 (note that the ECG energy along the two orthogonal dimensions of the decomposed space a„ and a., are proportional 
to the cori^sponding singular values ot--, or and <jj^. To g jaraniae that we deal with the morphologies rather than 
the energy differences, the decomposition space is rescaled to equalise the energies in bibth directions. 




50 



=[n, V?// yy^3 y^* 




Each Wj represents the reconstruction coefficients of the T wave df the j*, channel of the standard ECG 
[0063] The angle between different Wy vectors Is calcjiated; 
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5 g[0',180°] 

[0064] The smaller 6,y, the closer the reconstruction vectors of the i"^ and the j*'' ECG channels. It was observed in 
the data of this study that the T-wave morphology in V! is generally different than that of other channels, irrespective 

10 of any clinical background, mainly due tothe position oi the V1 electrode Figures 7a and 7b demonstrate the difference 
between a normal ECG and a HCM case. The projection of the T wave loop onto reconstruction vectors gives the 
particular T wave as cJaserved in the corresponding EGG lead. It is seen that reconstruction vectors for normal ECG 
are closely grouped (meaning sirriilar morphology), whereas they are dispersed (meaning different morphologies) for 
HCM patients Note that the reconstruction vector of VI In normal EGG is far from the others This is the reason why 

15 it is preferred to exclude the VI reconstruction vector from the calculations of the T wave IVIorphology Dispersion 
descrfptor (i e , the average of angles between all pairs of reconstruction vectors), to be explained in greater detail 
below 

[0065] The descriptor, Twave Morphology Dispersion fTMO), is defined, as. the mean of all .e,j excluding VI: 

20 

25 TMD IS a measure of the spatial T wave morphology variation 

[0066] Since the ascending and the descending parts of the T wave are luiown to correspond to different facets of 
the repolarisation process, descriptors TMDpre and TMOpg^^ which are defined in the same way as TMD with the 
'ascending' part of the T wave {tjg < f < fyp) used for T^fD^^ and the 'descending' part {irp^ti fj-^) used for T/MDj^f 
were also computed 

30 [0067] Both of the QRS and the T wave represented by S^pg, and S,- follow an approximate loop in the column 
space of u The orientation of the T wave loop Is determined by selecting the unit vector ej^,, with the maximum t 
wave energy, subsequently the unit vector e^a perpendicular to 67^-1 with the maximum energy and finally the unit 
vector ej3 perpendicular to both 67;^ andcj^- 

[0063] The descriptor Total Cosirie RJo_T(TCRT) is defined as the afvemge of the cosines of the angles between 
35 en and Sofls(/) (columns of Sq^^ for all i within [f.f^tf,^ This is a measure of the vector deviation between the 
depolarisation and the repolarisation waves 



TCRT = j;^cos{z(t,,,s^ (I))). 



[0069] Figures 8a and 8b show examples of the QRS and the T wave loop orientations in a normal ECQ and a HCM 
'>s case TORT measures the deviation between these two lobps It, in effect, measures the difference between the prop- 
agation directions of depolarisatirai and repolarisation waves Ivlegative TCRT values con-espond to large differences 
in the loop orientation in Figures 8a and 8b, it can b6 seen that the orientation of QRS loop and T wave loop are close 
to each other in the case of a normal patient, whereais they are far from each other in the case of an HCM (hypertrophic 
cardiomyopathy) patient 

so [0070] The T wave representation S j-, is de-normallsed, removing the effect of energy normalisation aiid its projection 
onto spa/i(eji, ej-g) is calculated orienting the T loop along e^^i The rectangular area encompassing the T loop is 
divided into n»:1 (in the present implementation n=4900, arbitrary constant) cells of equal size. The loop is closed 
with a straight line connecting the end points and re-samp:ed with equal sampling steps of the 2D space This re- 
sampling assures that there Is at least one sample in every cell that the loop passes through The numbers of cells in 

ss the loop area and in the outer area are counted and the descriptors. Percentage of the Loop Area (PL) and Percentage 
of the Outer Area (PO)afe calculated PL is the proportion between the number of cells inside the loop and the number 
of all cells Similarly, PO is defined for the cells outside the loop 
Note that PL + PO < 1 because there are cells occupied by the loop itself 
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[0071] The descriptor termed Lead Dispersion-1 [LD^], is caiculated using the decomposed, energy normalised and 
DC-compensated T wave S7-, as it was after tlie DC-compensation step. Sj^o (/) G spariu-^, Ug) V tj^< i< fj-^ follows 
a path m the (^-plane, the T loop The rectangular area containing that path is divided Into 100 equal cells (arbitr^ 
constant) LD^ is defined as the number of different cells that the path involves It serves as a measure of the temporal 

5 variation of the interlead reiations during T wave. 

[067i] The descriptor/ Z-earfOispere/on-^fLDg), IS defined similarly using th^ It isthenumber 

of different cells that the T loop passes through, excluding the straight line that was added to cbse the loop. The basic 
difference between LDi and LD2 is that LD, is calculated usingthe energy:nomialised decomposed signals (maximum 
EGG energy = 1 ), while LD2 is calculated using the original "decomposed signals.. 

10 [0073] These descriptors provide further metiiods of quantifying abnormalities of ventricular depolansation and re- 
polarlsation, and thus relate to further inventions disclosed herein 

[0074J The analytical system was implemented on a standard personal computer with Pentium 133IVIHz CPU and 
BOiVIB RAM, using Mallab Version 5.2.0 (The MathWorks Inc., 1998). The system was tested wjth standard 12-lead 
ECGs recorded by the MAC VU Electrocardiograph (Marquette Medical Systems, Milwaukee. Wisconsin, .USA) 10 

15 second recordings with 500Hz recording rate were acquired and the so-called median beat was obtained for each 
channel of the recording These median beats, sampled at 250HZ, were used in the analysts 
|0075] The inputs to system are 8, time-aligned median beats, each one being a representative of the EGG mor- 
. phology in the corresponding standard channel The main time consuming computations are the area calculations 
which involve a recursive algorithm, also implemented in Matlab. On the average, the compulation of all parameters 

20 forasingle recording takes 177 seconds If excluding the area related parameter&(PL and PO), the analysis takes 30 
seconds per recording Matlab commercially available library without any modification is presently used A purpose 
built library would increase the performance of the system considerably 
[007©] Three sets of ECG recordings were used in the study; 

2S (a) Standard resting 12-lead ECGs recorded in 11 00 normal heafthy subjects, 913 male, aged 3at1 2 years, range 

10-81 years 

(b) 10 supine resting ECGs were recorded In each of 76 normal healthy subjects, 37 male, aged 38 HO years, 
range 13-59 years In each individual, the serial ECGs were recorded one immediately after another using the 
same electrode attachments and without the subject moving during the whole recording session Data acquisition 

30 of each recording lasted 10 seconds and, including the slectrocardiograph handling, each series of the 10 ECGs 

was obtained wilhin 3 minutes 

(c) Using the same recording strategy, 1 0 supine resting ECGs were recorded in each of 63 patients with hyper- 
thropic cardiomyopathy (HCM), 44 male, aged 39* 14 years, range 12-71 years 

35 [0077] Using a research version of the commercial QT Guard software package (Marquette Medical Systems, Mil- 
waukee. Vifisconsin, USA), several conventional descriptors of repoJarisation patterns were calculated for each ECG 
for oomparison. This software aligns all beats with respebt to the Q wave imset and was programnned to use the 
dovmslc^e inflex tangent method to detect the T wave offset The following conventional QT inten/al and T wave 
parameters were considered: 

40 

i Global Qt Dispersion (G-QTd) = MaK (QT Intewal in 12 leads) - Min (QT inten^al ih 12 leads) 

ii Precordial QT Dispersion (P-QTd) = Max (QT Inten/al in 6 precordial leads) - Min (QT interval in 6 precordial 
leads) 

ill Area QT Dispersion (A-QTd) = All of the 12 leads are assumed to have the same T wave onset and offset 
45 points. The areas under the T waves are calculated and the points at which they reach 90% of the correspcxiding 
total area are marked for each lead The dispersion, (maximum minimum) of these markers over 1:2 leads is cal- 
culated 

iv Global J to Tpeak Dispersion (G-JTpd) = Max (J to Tpeak interval in 12 leads) - Min (J to Tpeak interval in 12 

leads) 

so V Corrected QT Intervat (QTc interval). Bazetl formula corrected maximum QT interval in all 12 leads 

Principal Component /Vialysis of 12 leadT waves is also incorporated in the QT Guard packages The B components 
with associated eigenvalues- are obtained. Each eigenvalue is a measure, of the significance of the correspionding 
component If s, denotes the eigenvalue associated with the i<h prineipai component, the foltowing descriptors are 

55 calculated. 
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PCA ratio 1 (PGA,)- . J xlQQ. 



PCA ratio 2 (PCAg) = ^ X 1 QO 



PCA ratio 3 (PCA3) = X 100 
Note that these pafameters. can also be calculated by SVD, assigning s,= a,-. 

[007S] To find out whether the new methods assess EGG qualities additiorial to.the. conventional parameters, the 
correspondence between the new and the conventlonai parameters was investigated. All new and conventional pa- 
rameters of 1100 ECGs acquired frorn normal subjects were used for this investlgatfon and to calculate the Pearson 
Product-^foment coixelation coefficient (Statistica Package, Release 5 1) between the new and the conventjonai pa- 
rameters and the ages of the subjects. 

[0079] The reproducibility of all the parameters was assessed based on the variation of the measurements between 
serlaf EGG recordings from the same individual of the populations of 76 normal and 63:HCM subjects The ratio of the 
individual range to the total range was calGulated for each patient and each parameter More precisely, for a fictional 
parameter X, this ratio R^^j tor normal subject j Is equal to 



inax(X;)- niinfA';) 
maxlxf )- mm{Xf } 



where N=76 for our data set of 76 normal healthy subjects and is the value of descriptor A'ln m'^ EGG of the subject n 
[0080] The values ^ - were obtained In a similar way and for each descriptor X, the means and standard deviations 
of values fl^^^ and were calcuiated and used to compare the reproducibility of all ttie descriptors 
[0081] We™fso investigated the univariate and multivariate distinction between normals and HCM, using all of the 
new and conventional descriptors The parameters were compared on the basis of the signiftcancis In discriminating 
these two groups and in terms of specificity and sensitivity. The populations of 76 normal aand 63 HGM subjects were 
used and for each dsscfiptor the mean values of 10 recordings were considered 

[0082] Individual parameters were firstly Used in a univariate analysis The normal and HCM groups were compared 
using non-parametric Manri-Whttney test implemented usinig an in -house written software according to 1he original 
description (sea Maim H B. and Whibiey D R (1 947). "On a test of whether one ortwo random variables is stochastically 
larger than the other", Ann Math Statistics, 18, pp 50-60 for more information in this regard) P-value < 0 05 was 

considered as statistically significant 

[0083] The Receiver^ Operator Characteristic (ROG) cun/es VMhich showthe dependency of specificity on sensitivity were 
calculated for each individual parameter using an in^aouse software pacl<age (see Hnatkova K- , Poioniecki J D: Gamm A. 
J MalikM (1994). "Computation of muliifaetorial receiver operator and predictive accuracy charactenstics", Comp Meth 
Prog Biomed , 42, pp. 147-156, for more infomnation in this regard). The area under the ROC cun/e (reported as a per- 
centage) was used to characterise the predictive power of each parameter Independ^t of fixed sensitivity levels. 
[0084] Multiple regression analysis was used to assess the relative performance of individual descriptors in discrim- 
inating the HCM from normal subjects. The dichotomi' limit of each parameter was set to the mean of the average 
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values of the normal and the HCM groups Multiple regression models of different orders were calculated (Statistica 
package, Release 5 1) in a backward stepwise manner, by excluding the least significant variable at each step until 
the P'Values of ali sun/iving parameters were below 0 05 

[00853 The most significant parameters, which were identified by the multiple regression analysis, were used for the 
calcuiation of multivariate f=10C curves. The combinations of two and four parameters were used with the decisbn 
rules of at least 1 of 2 positive and at least 2 of 4 positive, respective^- : 
[0086] The results obtained using the new and conven-lonal descriptors were as follows: 

[0087] Table 1 gives the Pearson Praduct-Moment correlation coefficients between the conventional and the new 
descriptors and the age of the patients. None of the new or conventional descriptors had a significant correlation with 
the age ( I r [ < 0 16 for all parameters): the absolutevaluaof ail the correlation coefficients between new descriptors 
were < 0.5 except for: TMD^MDp,,^,: 0.91 , TMD/TMDp^: 0.93, TMDp^/TMDp^: 0 79, PLyPO. -0 94. PULD^: -0 54, 

PO/LDg. 0 50 . . ■ . , 

[0088] The absolute values of the correlation coaffleients between the: conventional and the new parameters, were 
all < 0 5. except: TMD/PCAg: 0 552. LDa/PCAa. -0.562 

[0089] The reprdducibllity of the conventional descriptors was generally poorer than that of the new ones (Figure 9) 
with the exception of PCA^ and PCA2, that had reproducibilities similar to tiose of the new descriptors 
[0090] Table 2 shows the meEin values of ej^observed in 1 1 00 normal subjects and confirms the reasons tor excluding 
VI lead f rom spatial variation descriptors An increased spatial variation oi T wave morphology in HCM subjects was 
observed in all of the three spatial variation descriptors, i.e TMD; TMOp,eandTMDp<^,. The mean TCRT was negative 
for HCM subjects and positive for normals The PL was larger for normal subjects Whereas the PO was smaller The 
mean value of LDi (as well as of LD2) m nonnal and HCM groups were close to each other 

[0091] Table 3 shows the comparisons between the descriptors in nornnal and. HCM subjects While all of the de- 
scriptors strongly ditlereritiate between both groups, sonne descriplors haye substantially lower p^lues than others 
The .QTc inten/al. and P-QTd offer tiie most significant, uniwriate differfinliation among the conventional descriptors 
However, TGRT, TMD arid TMbposjOutperlormed all pf the coiiventional descriptors, while TMDp„ had a p-yalue close 
to that of QTc iriterval which is the best among the COTventional descriptors. 

[g092] Table 3 also Shows the area under the univariate ROC curve for each descriptor. The results confirm the 
statistical comparisons. TCRT TMD and TMOp^s, have areas above 90%, QTc interval has the largest area (S5 6%) 
among the conventional descriptors 

[0093] Table 4 shows the p-values oJ each parameter in a succession of multiple regression models of different 
orders The descriptors TCRT TMOp^, P-Qtd and QTc interval survived throughout the successive multiple regression' 
models perforrned in a backward stepwise- fash ion When excluding TMDp„ and TMDpost. 'he fihat significant param- 
eters surviving the backward stepwise multiple regression analysis were TCRT, TMD, P-QTd and QTc interval, with p- 
values of 5.75X10-8, 0 011, 5 93x10^, 8.61X10^6^ respectively In both cases, TCRT outperformed all of the other 
parameters in all orders of multiple regression analysis 

Tablet , 

Pearson Produci-Moment correlation coefficlsnts between the new descriptors, the conventional descriptors and 
the age of Ihe subjects 



CORRELATION COEFFICIENTS 





TMD 


TWIDp^t 


TMDp^e 


TCRT 


PL 


PO 


LD, 


LDs 


TMD 


1 GO 
















TIWDpos, 


0 91 


1 00 














TIWDp,, 


0 93 


0 79 


■ 1 CND 












TCRT 


-0 0,1 


0 04 


0 05 


1 00 










PL 


-010 


-0 16 


-0 18 


-0 08 


1 00 








PO 


0 09 


0 14 


017 


0 08 


-0 94 


1 00 






LDi 


0 08 


0 05 


0 07 


0 00 


0 14 


-017 


1 00 




LD2. 


.0 .30 


-Q.1.8 


r0 17 


.020 


"0 54 


0 50 


-0 12 


TOO 


C3-QTd 


0 08 


0 10 


0 01 


-0 03 


O.04 


-0 02 


-0 01 


-0 08 


P-QTd 


0 16 


0 16 


0 09 


000 


0 03 


-0 02 


0 00 


-0 15 


A-QTd 


0 10 


0 11 


0.00 


-0 08 


004 


-003 


001 


-0 13 
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Table 1 (continued) 



CORRELATION COEFFICIENTS 




TMD 


TWIDpos, 


TMDpte 


TCRT 


PL 


PO 




LDs 


C3-JTpd 


0.23 


0.19 


0.09 


-0.12 


0.23 


-0.22 


0.02 


-0.48 




0.25 . 


0.17 


0.22 


0.09 


0.06 


-0.05 


0 15 




PCAa 


0.55 


0.46 


0 46 


-0.14 


0.04 


. -0,06. 


0,05 


: -0.56 


PCAs 


0.16 


0 15 


014 


0.03 


-0.11 


0.1G 


0.11 


-0.12 


QTc 
interval 


0.07 


0 07 


0 08 


-0 05 


-0 17 


0 16 


0 04 


0 04 


AGE 


-0 06 


-0,11 


0.00 


.0.09 


. -0.06 


0 01 


■ 0,05 


0;07 
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AVERAGt e,j VALUES 



Channel 


J 


II 


VI 


Y2 


V3 


V4 


vs 


V6 


I 


0 000" 


6:799"±- 
79970 




22 47°± 
24.82° 


I0 65°± 
10.79° 


4.648°± 


4 451.°± 
4 603° 


5.9()()''± 
6.218 


II 




0 mr 




27 44"! 
25.79" 




7575"t 
9.010° 


4 70i°± 
6,317° 


4.180°+ 
6.015° 


VI 






0 000° 




m 








V2 








0(XH)" 


12 94'^± 
21.72^ 


20 7r± 

24 39" 


24 86"± 
24 97" 


27.65"! 
25 42° 


V3 










0 000" 


a 6K5"± 
K93r^ 


1 2 84°± 
10.70" 


I5 64'-'± 
12,31° 


V4 














4 406°± 
5,460" 


7 209°± 
8.129° 


vs 














0000° 


2,939°± 
3.893° 


V6 
















0,000° 



Tabic 2 Average 6,j vaiues in the set of 1 100 normal subjects The shaded boxes contain the 6i,vi and 
6vi , values that are greater than the others 
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Table 3 



The maan values and the standard deviations of all parameters, the Mann-Whitney Test results and the area under 

univariate ROC curves of the separation between normal and HCM subjects 


SEPARATION BETWEEN NORMAL AND ABNORMAL ECGs 


Parameter 


Normal n=76 


HCM n=63 


Mann-Whitney Test 


Area under ROC 
Curves 


Mean±SD 


Mean + SD 


P-\felue 


TMD 


10 72 + 4784 


41 10 + 26 85 


2 818X10-18 


90 1% 


TiyiDp,,,^ 


6 141 ±4462 


36 68 ±27 49 


2 289X10"''9 


91 1% 


TIWDp„ 


8 682 + 4585 


4214± 32 62 


1 605X10 13 


851% 


TCRT 


0 522 ± 0 274 


-0 351 ±0 522 


3 548X10-19 


90 9% 
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Tables (continued) 



SEPARATION BETWEEN NORMAL AND ABNORMAL ECGs 


Parameter 


Normal n=76 


HCM n=63 


Mann-Whitney Test 


Area under ROC 
Curves 


Mean + SD 


Mean±SD 


P-Value 




0 671 + 0.CD85 


0.608 ±0,142 


5.935X10-3 


64.3% 




0 273 + 0^072 


0 328 ±0.115 


3.051X103 


65:2% 




36 40 ± 1 ,163 


34-81 ± 3.157 


2.522X10-S 


71.8% 


— — 

LD2 


724 5 + 346 1 


604. 9 ±458.1 


6.787X10-'' 


67 4% 




19 97 + 11 .62 


36 55 + 18 85 


6.989X10-3 


77.5% 




1 0 79 ± 8 776 


27 87+18 69 


6 611X10-11 


80.6% 


A-QTd 


13 70±S.564 


24.38+12 23 


2127XlG-« 


76 8% 


G-JTpd 


32 53+12.18 


45 96 + 20,61 


2 463X10 = 


70 8% . 


PCAi 


680 0 ± 226 3 


481 4 + 245 8 


6 698X10-8 


76 7% 


PCA2 


15 56 + 5 162 


23 56+10 85 


9.886x10-7 


74 4% 


PCA3 


4.826 ±2, 373 


7 765 + 4.235 


6.603X1 0-s 


.78,4% ■ 


QTc Interval 


404.4 ±15:27 


435.1 ±25 50 


4i22X1.0ri4 : 


85.6% , 
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P- VALUES IN MULTIPLE REGRESSION ANALYSIS 





P - VALUES i n M u 1 t i p 1 e 
Regression Models 


Model Order 


All (16) 




10 


9 


8 




4 


PARAMETERS 


TMD 


0 769 














TMD^i 


0 272 


0 140 


0 174 






TMD„„ 


0 037 


0.00« 


0.004 


0 008 


0 001 


TCRT 


! 45x10' 


3,43x10^ 


5 39x10 ' 


1 07x10' 


2 24x10 " 


PL 


0 880 










PO 


0.563 


0 060 


0,032 


0 054 




LD, 


0 378 










LD2 


0 758 










G-OTd 


0.335 


0.065 


0,1 10 


0 156 




P.QTd 


0.016 


0.006 


0.006 


0.015 


7 07x10"' 


A-OTd 


0.558 










G-JTpd 


0.092 


0.074 


0 121 


0 166 




PCA, 


0.628 










PCAj 


0 397 


0 122 


0 062 


0 082 




PCAj 


0352 


0205 








3Tc interval 


6 06x10* 


5 05x10 ' 


194x10* 


1 93x10' 


8 57x10-* 



Tabic 4 Significance ievelis of ihfe parameters m different orders of multiple regression models, 
calculated in a backward stepwise fashion by excluding the least significant parameter at each 
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[0094] The muttivariate ROC curve involving the descriptors TCRT, TMOprg, P-QTd and QTc interval, had an area 
of 98 4%. Using the conventional and nawdescriptctfs separately in bi-variate ROC curves, we obtained areas of 95 5% 
for TCRT and IMD^,^. and of 91 6% for P-QTd and QTc interval 

[0095] Hence, TCRT, TMDp„. P-QTd and QTc interval are mutually rndependent separators of nonnal and HCM 
EGGs of which TCRT is by far the strongest. 

[0696] The neW paramejers proposed in this specification are defined using the decomposition space and aimed at 
the description of the temporal and spatial variations of ventricular repblarisation. The descriptors TMD, TMOp^^ and 
TIWDpos, reflect the interleatfmbrphblogical variations of the T wave patterns, that is the spatiar variations. The area 
related descriptors, that is PL; PC, LD^ and LDg, characterise the temporal variations. TCRT introduces the concept 
ofxomparing the global wavefront directions of the depolarlsation and repolarisation processes 
[0097] The original hypothesis was that, compared to normal ECGs, the spatial and temporal variation of T wave 
morphology are increased and the depolarisation and the repolarisation vectors are more different in pathological 
recordings, such as In HGM patients. The statistical comparisons of this study verify this hypothesis, 
[0098] The mean of PL is higher and PO is lower in normal than in HCM subjects This suggests that the T loop is 
relatively smooth and connected (not crossing Itself) in normal ECGs than .in HCM ECGs On the other hand, despite 
the significant difference, the lead dispersion parameters (LDi and LDg) have similar mean values for both the normal 
and HCM subjects This suggests that the loop itsialf Is not discriminative The loop lengths were similar in both groups 
Th e discrimination by PL and PO seams to be due. to the disconnected and narrow loop (inner area is similai- to a strip), 
rather than due to an increased irregularity of its shape 

[0099] The change of sign of TCRT between the normal ahd the HCM subjects provides a clear distinction between 
the two groups (the negativity of TCFIT shows an increased deviation) This is in agreement with the original hypothesis 
The repolarisation and depolarisation waives do differ in tenns of their principal direction in a 3-dimensional time-or- 
thogonal space Since the mean difference between normal and abnom;tal ECGs ts-52° , the descriptor does riot merely 
reflect T wave inversion that would result |n thedffferencs nearto leO'' - 

[0100] TCRT, TMDpre, P-QTd and QTc interval were the only parameters that survived throughout the backward 
stepwise multiple regression analysis comparing normals and HCM subjects and TORT was the best throughout the 
test This verifies that the new spatial variation parameters and TCRT are very potent descriptors of repolarisation 
abnonnalities 

[0101] In this preferi-ed embodiment, all of the new descriptors are defined using the decomposition space This 
provides an inhcrcnl immunity to noise and avoids ttie inaccuracies associated with time domain measurements, that 
are common in QT interval related descriptors that depend on T wave offset determination The independence of time 
domain measurements makes the new descriptors highly reproducible, which is veiy important for their potential clinical 
applicability Among the conventional parameters, only PCA, and PCAg have a reproducibility in the same order which 
is again due lo avoiding the time domain measurements 

[01 02] The weak correlation between the new and the conventional parameters shows that the new concepts quantify 
different properties of the ventricular repolarisation Furtiiermore, the two new concepts, thatis flie spatial and ternporal 
variation can be identified by strong correlations within each group and by weak correlations across the groups The 
concept of TCRT is different from both spatial and temporal variations and the descriptor does not correlate strongly 

with any other descriptor 

[0103] The relatively poorer reproducibility of the T loop area related parameters is due to the algorithmic problems 
An open loop may result from baseline wander, as well as ST-segment elevation/depression A sliaighl lino was used 
to connect the ends of the loop, which is not necessarily the best approach An alternative may be to connect the ends 
of the loop and Its centre of gravity or to transform the w, Uj-plane creating a closed loop. It \s, also possible that the 
loop crosses itself, resulting in more than one surrounded area in the study the inner area was defined as Hie closed 
area neighbouring the beginning of the loop, ignoring the "pockets" It is an open question vyhether the existence and/ 
or the area of ttiese 'pockets' is of any significance. The poor perfomianceof the loop related descriptors In differen- 
tiating normal and ^bnonnai ECGs may well be due to these problems 

[0104] The arbitrary choice Of constants. Used in dividing the plane of ECS into equal size cells, have an influence 
on LDi , LD2, PL and PC calculations as well as in approximate T wave offset detection They define the precision of 
these descriptors. Increasing these constants would increase the precision at the cost of Increased computation time 
However, the precision is also restncted by the ECG sampling rate which determines the smallest distance between 
two consecutive ECG vectors Unreasonably decreased cell size [increased constants) would also degrade the per- 
formance of T wave offset detection 

[0105] On the other hand, the T wave onset/offset definitions may have an influence on the temporal variation de- 
scriptors but do not affect the others Setting the constant p. in T wave offset detection to 3 is an appropriate choice 
The algorithm readjusted that value in 91 of 1100 nonnal ECGs The QRS onset/offset definitions , on the other hand, 
are robust and able to handle wide QRS complexes However, the choice of 70% threshold in determining the region 
of QRS used in TCRT calculatkjn is important A too low threshold may result in a too general estimation of the QRS 
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loop orientation, whereas a high tlireshold may misinterpret the orientation of tlie depolarisation wavefront vector. 
[010S] Only the principal direction of the EGG vector during T wave, Bj-, was used in TCRT Galculation, thal is the 
other lower energy components, Cy-.s and 67-3 were ignored The average ratio of the energy along the second ccm- 
ponent to that of the first was 0, 1 4for normals and 0.22ior HCM patients. This shows that the T toop generaliy resembles 
s a narrow ellipsoid and it is the direction of this loop that is of interest. Using 67-3 and/or Ci-^g would not improve this 
concept of TCRT, mainly due to a decreased rioise inimunity. There is no annbiguity in the Ot-:, definition because the 
DC-compensation ensures that ©7^1 Kas the: correct sign. ^ .V 

[010^ In conclusion, therefore, it can be seen that the new descriptors of repolarisation patterns descrtbed. in this 
specification have seVeralimportant qualities. 

10 : 

• All of the new descriptors can be assessed in a minimum dimensional space constructed by SVD of 12-lead ECG 
This provides a buitt-in immunity to npise. 

• None of the new descriptors require accurate time domain interval measurements This makes them more repro- 
ducible than the conventional QT ihten/al related descriptors 

)S • The new descriptors assess different ECG qualities than the conventional parameters 

• The spatial variation and wavefront direction descriptors can discriminate between normal and abnormal EGGs 
- - substantially better than the conventional descriptors TTie wavefront direction descnptor (TCRT) is by far the 
/ : strongest of all considered in this study 

zo [0108] In a second case study, the effect of changes In the autonomic tone of a patient was studied using known 
descriptors and some of the new descriptors described in the first case study The findings were as follows. . 
[0109] In the early thirties, Wilson efal (Wilson FN . MacleodAQ, Barker RS., Johnston ED , The defermiriation 
and the significance of the areas of the ventricular .deflections of the eleclrocardipgram A. Heart J 1 934; 10/46-61 , 
Wilson F N , Macleod A G , Baker RS.^ The T Deflectbn of the Electrocardogram. TR A Am Physicians 1931; .46.29) 

zs ihtroduoed the concept of algebraic sum of the areas under the ventricular deflections of the aiectrpcardiogranri <the 
net QRST area) The resulting vector quantity, called 'the ventricular gradient" (VG) was proposed as independent of 
the sequence of ventricular activation as long as the ventricular recovery properties remaned constant Therefore it 
was believed that the VG could help distinguish T wave changes following changes in the activation pattern ('secondary' 
T wave changes) from those due to myocardial damage ('primary' T wave changes) 

30 [0110] The original idea did not fivolve inlo a clinically useful tool both becau.'se of technical difficuitles with the meas- 
urement, and because of data challenging the quantitative independence of the activation sequence, As discussed in 
relation to the first case study, a new descriptor of the wavefront direction of repolarisation has been proposed It 
quantifies the differaice in. the global direction oflhe wavefronts of.thadepolarisation and of the repolarisation as an 
average of the cosines of angles between mam depolarisatipn and repolarisation vectors in a minimum djmensbna] 

35 subspacederivedfromtheindependentleadsof the electrooardtogram(t(italcosineRtoTTCRT) Although measured 
in the optimised 3D space that contains the maximum energy of al| ECG leads, this new descriptor advances the 
classical concept of VG 

[0111] Ventricular gradient is a vector which gives the direction and magnitude of the electrical forces produced by 
a lack 0! uniformity jn the .duration oflhe excited state, it points from the region in which the average length of systole 

40 is greatest, toward the region in which it is least Once considered a 'fundamental quantity in electrocardiography' the 
VG was gradually forgotten due to both uncertainties about the validity of its concept as well as to technical difficulties 
with Its manual calc jiation from QRS-T time integrals before the personal computer era Today the VG is hardly men- 
tioned in modem textbooks of electrocardiology with comments such as that 'the most exciting thing about the ven- 
tricular gradient is its name' The TGRT has been demonstrated to be more reproducible and to separate normal from 

45 abnormal ECGs better than several repolarisation parameters including the dispersion of the QT-end and JT-peak 
intervals and the corrected QTc intenral TCRT has also been found to contain independent predictor value of mortality 
and arrhythmK complicalions after myocardial infarction. 

[tli12] These obsen/ations suggest thafthe concept: 0!f VG might have been neglected prematurely Having this in 
mind, several studies were initiated researching the properties of VG and of its modern andmore precise counterpart 

SO In the present study the effects of basic autonomic provocations on VG were investigated The study assessed the 
effect of postural changes and autonomic provocative manoeu^'res on the direction and magnitude of the spatial VG 
and on TGRT in healthy subjects We compared the effects of postural and autonomic provocations on VG descriptors 
with the effect on conventional repolarisation parameters, namely theOT interval duration and QT dispersion 
[0113] The study population consisted of 40 healthy subjects. 31 male, median age 33 years, mean age 33.1+7 3 

55 years, range 1 8-56 years, with no history of heart disease and with nomial resting 1 2-iead ECG None of the subjects 
wias takingany autonomically active medication and before the test, the subjects were instructed to refrain from smoking 
and from caffeine intake 
[0114] The procedure was as follows. 



17 



EP 1038 498 A2 



[0115] Following 10 minutes supine rest in a comforting tQmperatured and dimmed room with a bw level of back- 
ground noise, the subjects perfornned tlie following testS- 

[0116] Postural changes (32 subjects), resting supine position, followed by sitting, unsupported standing, supine, 
and standing position, 4 minutes in each position {total of 20 minutes), with abrupt transition between the separate 

s positions. 

[0117] Valsalva manoeuvre (30 subjects): forced e)q)iration into the mouthpiece of mercury manometer maintaining 
a constant pressure of 40 to 50 mm Hg for up to 1 minute. Each manoeuvre was perfomied 3 times in supine position 
with 4 min periods between the tests, and 3 times in unsupported standing p&sition also with 4 minutes periods between 

the manoeuvres. _ . 

70 [01 1 8] Sustained Handgrip (8 subjects); the maximum force of contraction was determined in each individual follow- 
ing which each subject maintained 30% of maximum force for 5 min. Each test was performed twice in supine position 
with 4 minute periods between the tests and twice in unsupported standing position again with 4 minutes of rest between 
the tests 

[01 1 9] The data were recorded using continuous 12-lead digital EGG (250 Hz sampling rate, 1 2 bit A/D conversion) 
is for the complete procedural section for each subject without any toss of signal using a digital recorder (SEER MC, 
Marquette Medical Systems, Milwaukee, Wisconsin, USA) EGG data were stored in separate 10 sec portions The 
individual tests were organised in a synchrony with the recorder in order to identify each 1 0 EGG sample within a 
specific phase of each test 

[0120] From each lead of each 10 sec EGG sample, the so-called 'median beat' was constructed representing the 
20 ideal QRST complex of the given EGG Data analysis Of EGG patterns was based on these median complexes 

[0121] The magnitude and angle of spatial VG were calculated in the ifollowing ways From the median EGG beat, 
the area of the QRS complex and of the T wave in each of the 12 leads were calculated using the EGG Research 
V^oritstation Software Package Version 1 0, by Marquette G E. (Milwaukee, Wl, USA). The QRS and T wave areas in 
the orthogonal X, Y, Z leads were derived from the areas in the 8 independent leads (I, II, V^, V2, V3, V4, Vg and Vg) 
55 using the inverse Dower matrix. The riiagnitude of the spatial VG was calculated as 

VG^ = ((QRS^ + t/ + (QRSy -1- t/ + (QRS^ 4- t/)'* 

30 where ORS^,, and T^, are the areas of the ORS complex and of the T wave in the orthogonal lead W, respectively The 
angle VGa of the spatial VG was calculated as the spatial angle between vectors originating in the centre of the co- 
ordinates with the final points of [QRS^.QRSy, QRSJ and[T^.Ty.TJ, respectively 

[0122] As described in detail in the descriptbn of the first case study, the total cosine R_to^T (TORT) was calculated 
The median beat was represented in a minimum dimensional subspace using singular value decomposition of the 
3S Standard 12-lead EGG The TCRT was defined as the average of the cosines of the angles between the main QRS 
and T vectors in the 3-dlmensiqnaJ reconstmeted subspace. In effect, tGRT measures the difference between the 
propagation directions of the depolarisatioh and repolarisatron waves, with smaller (and negative) values Gorrespondin g 
to greater difference between the two wavef fDnt directions. 

[01 23] The median beats of all ECGs were analysed automatically using the QT Guard software package (Marquette 
40 G.E ) A common onset of the Q-wave in all leads was identified and the offset of the T wave in each lead was determined 
by the downslope inflex tangent method For the purpose of this study the maximum QT inten/ai, the global QT dis- 
persion (QTd, maximum QT interval in 12 leads - minimum QT inten/al in 12 leads) and the RR inlerva! were taken 
from the results provided by the QT Guard package 

[01 24] The first 2 minutes of the supine rest recordings were ignored in order to achieve f u lly stabilised steady state 
45 The recordings obtained during the final 8 minutes of lhe 10 minutes ol supine rest were used to derive baseline values 
for each parameter and to investigate their mutual correlaticn as well as correlation with heart rate To investigate the 
correlation, averaged supiiie resting values of Individual subjects were corisidered. 

[01 25] The mean values of each parameter for the separate positions and autonomic manoeuvres were calculated 
and compared by paired t-test and one-way within subjects (repealed measures) analysis of variance (ANOVA) with 
so post hoc comparisons using Scheffe test (Statistica, Version 4.00) All values are expressed as mean ± standard error 
of the mean (SEM) Statistical significance was defined as p<0 05 

[01 26] The results of the study were as follows. The correlation coefficients between VGa, VG, TCRT and RR.interval 
during steady-state supine position are shown in Table 5. There was a significant correlation between the angle and 
the magnitude of the VG, and between the angle of the VG and TCRT While both VGa and VGm were significantly 
ss correlated to the RR interval, there was no significant correlation between TCRT and RR (Table 5) 

[0127] Postural changes significantly decreased both VG and TCRT. While VGm and TCRT were significantly de- 
creased In sitting and further in standing position position, VGa was increased in sitting and was further increased in. 
standing position (Table 6) 
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[0128] As expected, the RR interval was significantly shortened in sitting and further in standing position. The max- 
irrium QT interval followed the changes of the RR interval (see Figures 10a and 10b). 

[0129] VGa.VC3m and TCRT were abruptly changed with transition from supine to sitting and from sitting to standing 
position, Ther^e was a general tendency for TCRT to be changed more abruptly when assuming each new position, 

5 compared to both VGa and VGm. 

[0130] VGm and TCRT were significantly decreased and VGa was significantly increased during the strain phase of 
Valsalva manoeuvre compared to preceding resting period both in supine and in standing position 
[0131] The RR interval were significantly decreased during the strain phase of Valsalva both supine and standing. 
QT max was sligthly but statistically significantly shortened during Vfeilsalva in standing (364±4 vs 371 ±5 ms, p=0 02) 

to but not in supine 
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position (Table 6). 

[0132] QT dispersion was not changed significantly during Valsalva 

[0133] VGmand TCRT were slightly but significantly increased by handgrip in supine position compared to the pre- 
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ceding rsstjng period (57.8+4.1 vs 56 1±3 9 mVms, p<0.0001 for VGm and O.SatO.Oe vs 0.61±0.06, p=0.0007 for 
TCRT) Both descriptors were not changed significantly by handgrip in standing position VGa was slightly but statis- 
tically significantly reduced by handgrip In supine (40 0+3. 1 vs 40 8+3 1 , p=0 01 ) and increased in standing position 
compared to preceding resting period (52 8+3 8 vs 51 .8±3 7, p=0 027) 

[01 34] Neither the RR interval, nor QTd were changed significantly by handgrip compared to preceding resting penod 
in supine as well as in standing positbn. However, QTmax was slightly, but significantly shortened by handgrip.in 
standing position 370+5 vs 374+5, p=0 06). 

[01 35] The aim of this study was to examine the effect of established autonomie tests on the spatial ventricular 
gradient and on a new descriptor of the wave direction of depolarisation and repolarisation, TCRT. 
[0136] The main finding is that both the magnitude and the angle of the spatial veintrieular gradient, as well as the 
new wave direction descriptor, TCRT, react sensitively and rapidly to changes induced by postural aid Valsalva ma- 
noeuvres As expected from Weir mathem.atical relation (i.e the algebraic sum of two vectors decreases as the angle 
between them increases) VGm followed the opposite trend to VGa and was significantly reduced m sitting compared 
to supine and was further reduced in Standing position 

[0137] Even more impressive than the absolute magnitude was the speed of the change of VG and TCRT during 
postural changes and Valsalva manoeuvre Statistically significant changes in TCRT were detectable already in the 
first or second 10-sec median EGG beat after standing up from supine or sitting position, or lying down from standing 
position (Figures lOa and 10b) Considering the fact that during the first 10-sec recording in each new position, the 
actual postural change took place (taking approximately 3-4 seconds), it is apparent that TCRT is an extremely rapidly 
responding parameter TCRT appears to be able to respond more rapidly to autonomic modulations than both VGa 
and VGm 

[0138] Previous studies have described significant reduction of VG^: In standing position compared to supine, as 
well as duringthe strain phase of Valsalva manoeuvre. The effect of all phases of Valsalva manoeuvre on repolarisation 
descriptors; howeveir, c&n bis assessed precisely only on beat-to*eat basis analysis 

[0139} Itis difficult to comrhent on the differences in the effect of Valsalva manoeuvre and sustained handgrip on the 
wavefront direction descriptors Although handgrip in supine position statistically significantly increased VGm (56. 1±3, 9 
vs 57 8+4 1 mVms, p<0 0001) and TCRT (0 61±0 06 vs (1 63+0 006, p=0.0007) and decreased VGa (4Q.8+3.1 vs 
40 0±3 1 , p=0.01 ) the differences were much smaller than those induced by Valsalva manoeuvre.and postural changes 
The handgnp test, however, is known to be of limited sensitivity and specificity, 

[0140] We found statistically significant correlation (r = -0 78, p <.0 0001) between VGa and TCRT during steady- 
state supine conditions The spatial VGa and TCRT appear to quantify the same physiological phenomencMl, namely 
the difference in the spatial direction of the. wavefronts of the depolansation and the repolarisation, in a three-dimen- 
sional space, u J 
[0141] VGa was positively correlated with heart rate (r = -0 39, p=0,026 for the VGa/RFi relation), On the other hand, 
although TCRT and the heart period clearly followed the samelrehd during postural, changes and . Valsalva manoeuvre 
(see Figures 10 and 11 ), there was no correlation between the two parameters during steady state supihe position (r 
0 09 p = 0 61 ) This suggests that rather than JCHT being driven by the heart rate, both parameters are under the 
control of a common factor'^operating during autonomb perturbations but not during steady state resting conditions 
[01 42] Previous studies have found significant correlation between VG and the heart rate . On the basis of this relation 
VG found limited use as a sensor for rate-adaptive pacing The gradient Is increased by an increase in pacing rate, 
and IS decreased by exercise or other stress during fixed rate ventricular pacing. Thus, at least in theory, the gradient 
can be implemented in aclosed loop rale-adaptive feedback system, the Bradierit js decreased by physical or emotional 
stress which leads to increase of the pacing rate (at a controjiable speed) and the latter increases the gradient to the 
resting level. The gradient was measured by the time Integral (the area) under the evoked ft wave Iniportant limitations 
of the VG measured in this way prevented the widespread use of VG as a rate-responsive sensor: although it is a 
rapidly responding sensor, the maximum pacing rate is reached very early during exercise, so the proportionality of 
the ratefesponse is very nodoiats The ventricular gradient is also directly related to the ventricular mass and thick- 
ness, therefore changes m the ventricular geometry in upright posture may lead to paradoxical changes in heart rate 
[01 43] The data from this study suggest that the role of VG and especially of TORT as rate-responsive sensors in 
implantable antiarrhythmic devices, and possibly aiso for automatic detection of potentially arrhythmogenic autonomic 
modulationsare potential applications of these descriptors Toachievethisinacommercial product, it may be necessary 
to make certain modifications such as simplyfying the measurement algorithm of TCRT and to calculate it from rntrap- 
ardiac Isads instead of standard surface leads 

[0144] In some recent studies it was found that the beat-to-beat variability of the VG was significantly increased by 
myocardial ischemia (see, for example Horinaka S, Yamamoto H, Tabuchi T Takada M, Akabane T, Onoda M, Yagi S 
Ventricular gradient variability New ECG method for detection of ischemic heart disease J Electrocardiol 1995; 28. 
177-183 and Horinaka S,. Yamamoto H Enhancemenf of ventricular gradient variability during apute myocardial 
ischemia Int JCardlol 1998, 65.173-180) In another re=ent study it was found that an abnomial T wave axis in either 
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the frontal or horizontal plane was a strong and independent prsJictor. of fatal and noPTfatal cardiac events in the 
general population older than 55 years: TCRT is able to observe such an abnornial -T-wave axis. 
[0145] From the results of this study, it was found that the maximum QT interval followed the changes of heart rate 
during postural changes and was significantly decreased in sitting compared to supine position and was further de- 
s creased in standing position. Although statistically significant, the changes of several millisecond^ of QT max in standing 
position by Valsalva and handgrip are hardly of any clinical significance. 

[0146] Although QT dispersion was significantly decreased in sitting and further in standing position compared to 
supine, the values are largely overlapping. It has been shown in the past that QT dispersion was significantly increased 
in standing compared to supine position in patients with syndrome X, while in one study a significant eflect of posture 
^o on QTdispersion in healthy subjects was notfound The low reproducibility of both automatic and manual measurement 
of the QT dispersion is well documented 

[0147] The role of the autonomic neivous system in ventricular arrhythmogenesrs has recently been heavily inves- 
tigated (sec, for example. La Ftovere MI, Bigger Jr JT Marojs Fl, Mortara A, SchwarU PJ Baroreflex sensitivity and 
heart rate variability in predictioi of total cardiac mortality after myocardial infarction Lancet 1 988, 351 .478-484, Hohn- 

15 loser SH, Kllngenheben T, Loo A, Hablawetz E, Just H, Schwartz PJ . Reflex versus tonic vagal activity as a prognostic 
parameter In patients with sustained ventricular tachycardiaor ventricular fibrillation Qirculation 1994, 89.1068-1073, 
arsd Schmidt G, Malik M, Barthel P, Schneider R, Ulm K, Rolnitzky L, Gamm AJ, Bigger Jr JT, Sehomig A Hear-rate 
turbulence after ventricular prematura beats as a predictor of mortality after acute myocardial infarctlort Lancet 1999, 
353: 1 390- 1 396) The autonomic modulations of the dynamlcity of ventricular repolarisation have been shown to contain 

20 diagnostic and prognostic information (see, for example, Maison-Blanche P, Coumel P Changes in Repolarization 
Dynamlcity and the Assessment of the Arrhythmic Risk PACE 1 997, 20[Pt ll]..261 4-2624, NollO G, Speramza G, Grasso 
R, BonaminI R, Mangiardi L, Antolini R; Spontaneous besat-to-beat vari^ility of the ventricular repolarisation duration 
J Electrocatdiol 1992, 25.9-17; Sapdrone G, Torzillo D, Fundarb.C, Porta A, Danna P, Pplese A, Malliani.A, Lombardi 
F Spectral Analysis of RR and R-T Variabilities in Patients With Coronary Artery Disease. A N E. r998; .3(3);237-243; 

25 Gang Yi, Xiao-Hua Guo, Reardwi M, Gallagher MM, Hnatkova K, Camra AJ, Malik M Circadian Nferiation of the QT 
Interval In Patients With Sudden Cardiac Death After Myocardial Infarction Am J Cardiol 1 99B, 81 950-956; and Horns 
E, Marti V, Guindo J, LagunaP, VinolasX, Caminal P, Elosua R, Bayesde Luna A. Automatic Measurement o( corrected 
QT interval in Holier recordings. Comparison of its dynamic behaviour in patients after myocardial mffiirction with and 
without life-thraatening arrhythmias Am Heart J 1997,134.181-7 However, the clinical valiis of many approaches to 

30 the assessment of the repolarisation dynamlcity is unclear Reliable descriptors to dotoct and quantify the autonomic 
modulations of ventricular repolarisation, would he extremely useful. The sensitivity of VG and of the new decriptor 
TCRT for detecting autonomic modulations, could be translated into prognostic power for major cardiac events As 
mentioned above, the preliminary results indicate that such use of these descriptors is certainly possible 
[0148] the conclusion, from this study was that both VG and TCRT are superior to QT max and QTd in detecting 

3S autonomic modulations o1 myocardial repdiansation induciad by postural changes and Vklsalva manoeuvre The sen- 
sitivity of TCRT and the dynanriicity of its reaction to postural and autonomic modulations suggest that it would have 
particularly useful application for rate-adaptive sensors in implantable antiarrhythmic devices Given its measurement 
could be simplified and made possible from itracardiac leads it could be used for automatic detection of autonomic 
modulation of ventricular repolarisation or other potentially anythmogenic factors in implantable devices 

40 [0149] There now follows a description of preferred implantable devices, such as a pacemaker (see Figures 12 to 
14) and a pacemaker cardioverter-defibrillator (see Figures 15 and 16), which could use TCRT to improve pacing to 
act as a monitor for detecting ventricular repolarisation abnormalilies or even act as an alarm to vrarn of autonomic 
conditions that pose a risk to the patient 

[0150] Figure 12 is a isimplified schematic view of one embodiment of an implantable medical device ("IMD") 10 of 
45 the present invention I MD 1 0 shown In Figure 1 2 is a pacemaker comprising at least one ol pacing and sensing leads 
1 6 and 1 8 attacfied to hermetically Sealed enclosure 1 4 and implanted near human or mammalian heart 8 Pacing and 
sensing leads 16 and 18 sense electrical signals attendant to ffie depolarization and re-polarization of the heart 8, and 
further provide pacing pulses foi- causing depolarizatkjn of cardiac tissue in the vicinity of the distal ends thereof Leads 
16 and 18 may have unipolar or bipolar electrodes disposed thereon, as is well known in the art Exarnples of IMD 10 
so include implantable cardiac pacemakers disclosed in U 8 Patent No. 5,158,078 to Bennett et al , U S Patent No. 
5,312,453 to Shellon ei al or U S Patent No. 5,144,949 to Olson, all hereby incorporated by reference herein, each 
in its respective entirety 

[0151] Figure 13showe connector headermodule12and hermetically scaledenclosure 14 of IMD 10 located in and 
near human ormammalian heart 8 Atrial and ventricular pacing leads 1 6 and 1 8 extend from connector header module 
£5 1 2 to the right atrium and ventricle, respectively, of heart 8 Atrial electrodes 20 and 21 disposed at the distal end of 
atrial pacing lead 16 are located in the right atrium Ventricular electrodes 28 and 29 at the distal end of ventricular 
pacing lead 18 are located In the right ventricle 

[0152] Figure 14 Shows a block diagram lllustratihg the constituent components of IMD 10 in accordance with one 
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embodiment of the prGsent invention, where IMD 10 is pacennaker having a microprocessor-based architecture. IMD 
1 0 is shown as including activity sensor oraccGlerorneter 11 . which IS preferably a piezoceramic.accelerometerbpn 
to a hybrid circuit located inside enclosure 1 4 Activity sensor 11 typically (although not necessarily) provides a sensor 
output that varies as a function of a measured parameter relating to a patient's metabolic requirements Additional 
information based on the patient's autonohiic tone could be provided here to improve the pacing. For the sake of 
convenience, IMD 10 in Figure 14 is shown with lead 18 only connecied thereto; similar circuitry and.connections not 
explicitly shown in Figure 14 apply to lead 16. ^ /^ -^^..^ 

101531 IMD 10 in Figure 14 is most preferably prograrTirnable oy means of an. external programming unit (not. shovm 
in the Figures) One such programmer, is the commercialiy available Medtronic Model 9790 programmer, which is 
microprocessor-based and provides a series of encoded signals to IMD 10,. typically, through a progrannming head 
which transmits or telemeters radio-frequency (RF) encodedsi9nals to IMD 10. Such a telemetry system is described 
in U S Patent No 5,312.453to Wybemy etal , hereby incorporated by reference herein in its entirety. The programming 
methodology disclosed in \%borny et al.'s '453 patent is Identified herein for illustrative purposes only Any of a number 
of:suitable programming andlelemetry methodologies known in the art may be employed so long as the desired infor- 
mation is transmitted to and from the pacemaker ^ , , 
[0154] As shown in Figure 14. lead 18 is coupled to node 50 in iMD 10 through input capacitor 52 Activity sensor 
.or accelerometer 1 1 is most preferably attached to a hybnd circuit located inside hermetically sealed enclosure 14 of 
.eiMD 10 The output signal provided by activity sensor 11 is coupled to input/output circuit 54. Input/output circuit 54 
;^contains analog circuits for interfacing to heart 8, activity sensor 11, antenna 56 and circuits for the application' of 
stimulating pulses to heart 8 The rate of heart 8 is controlled by software-implemented algorithms.stored m microcom- 
! puter circuit 58 

rOISSI Microcomputer circuit 58 preferably comprises on-board circuit 60 and off-board circuit 62 Circuit 58 may 
.correspond to a microcomputer circuit disclosed in U S Patent No. 5,312.453 to Shelton at al , hereby incorporated 
by reference herein- in its entirety On^faoard circuit 60 preferably includes nnicroprocessor 64, system clock circuit 66 
and on-boart! RAM 68 and ROM 70 Off-board circuit 62 preferably comprises a RAM/ROM unit On board circuit 60 
and pff-board circuit 62 are each coupled by data communication bus 72 to digital control lerAimer circuit 74 Micro- 
computer circuit 58 may comprise a custom integrated circuit dovice augmented by standard RAM/ROM .Gomponents 
[0156] Electrical components shown in Figure 1 4 are powered by an appropriate implantable battery power source 
76 in accordance with common practice in the art For the sake of clarity, the coupling of battery power to the various 
components of I MO 1 0 is not shown in the Figures Antenna 56 Is connected to input/output circuit 54 to permit uplink/ 
downlink telemetry through RF transmitter and receiver telemetry unit 78. By way of exarripte, telemetry unit 78 may 
correspondtothat disclosed in U S PatentNo, 4,566.063 isisuedtoThompson etal .hereby incorporated by reference 
herein in its entirety, orto that disclosed in the above-referenced '453 patent to Wyborny etal It is generally preferred 
that the particular programming and teleirietry scheme selected permit the entry and storage of cardiac rate-response 
parameters The specific embodiments of antenna 56, input/output circuit 54 and telemetry unit 78 presented herein 
are shown for illustrative purposes only, and are not intended to limit the scope of the present invention 
[0157] Continuing to refer to Figure 14, Vref and Bias circuit 82 most preferably generates stable voltage reference 
and bias currentsforanalogcircuits included in input/output circuit 54 Analog-to-digltalconverter (ADC) andmultiplexer 
unit 84 digitizes analog signals and voltages to provide "real-time" telemetry Intracardiac signals and battery . end-of- 
life (EOL) replacement functions. Operating commands for controlling the timing of IMD 10 are coupled by data bus 
72 to digital controller/timer circuit 74. where digital timers and counters establish the ovarall escape interval of the 
IMD 1 0 as wall as varbus refractory, blanking and other timing windows for conlrolling the operation of peripheral 
components disposed within input/output circuit 54 

[0158] Digital controller/timer circuit 74 is preferably coupled to sensing circuitry, including sense amplifier 88. peak 
sense and threshoW. measurement unit 90 and comparator/threshoW detector 92 Circuit 74 is further preferably cou- 
pled to electrogram (EGM) amplifier 94 for receiving amplifiedand processed signals sensed by lead 1 8 Sense amplifier 
88 amplifies sensed electrical cardiac signals and provides an amplified. signal to peak sense and threshold measure- 
ment circuitry 90, which in turn provides an indication of peak sensed vqllagesandmeasured sense amplifier threshold 
voltages on multiple conductor signal path 67 to digital cQntroller/limer circuit 74. An amplified sense amplifier signal 
is then provided to comparator/threshold detector 92 By way of example, sense amplifier 88 may correspond to that 
disclosed in U S Patent No 4,379,459 to Stein, hereby incorporated by reference herein in its entirety 
[0159] The electrogram signal provided by EGM amplifier 94 is employed wrfien IMD 10 is being interrogated by an 
extemal programmer to transmit a . representation df a cardiac analog electrogram See, for example, U S Patent No 
4 556 063 to Thompson St al. . hereby incorporated by reference herein in its entirety Output pulse generator 96 pro- 
vides pacing stimuli to patient's heart 8 through coupi^ig capacitor 98 in response to a pacing trigger signal provided 
by digital controller/timer circuit 74 each time the escape inteival times out, an externally transmitted pacing command 
is received or in response to other stored commands as is well known in the pacing art By vvay qf example, output 
amplifier 96 may correspond generally to an output amplifier disclosed in U S Patent No. 4,476,868 to Thompson. 
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hereby incorporated by reference herein in its entirety. 

[01 60] The specific embcxJimsnts of inpirt amplifier 88, output amplifier {output pulse generator) 96 and EGM amplifier 
94 identified herein are presented for illustrative purposes only, and are not Intended to be limiting In respect of the 
scope of the present invention. The specific embodiments of such circuits may not be critical to practicing some em- 
bodiments of the present inventicn so long as they provide means for generating a stimulating pulse and are capable 
of providing signals indicative of natural or stimulated contractions of heart 8. 

[0161] In some preferred embodiments of the present invention, IMD 10 may operate in various non-rate-responsive 
modes, inciuding, but not limited to, DDD, DDL WI. VOO and WT modes. In other preferred embodiments of the 
present invention, IIVID 10 may operate in various rate-responsive, including, but not limited to, ODOR, DDIR, VVIR, 
VOOR and VVTR modes. Sonrie embodiments of the present invention are capable of operating in both non-rate- 
responsive and rate responsive modes. Moreover, in various embodiments of the present invention IMD 10 may be 
jsrogrammably configured to operate so that it varies the rate at which it delivers stimulating pulses to heart 8 only in 
response to one or more selected sensor outputs being generated One such response could be generated through 
measurement of TCRT. This could be in response to changes in the autonomic tone of the patient or changes in direction 
of the repolarisation wavelront Numerous pacemaker features and functions not explicitly mentioned herein may be 
incorporated into IMD 10 while remaining within the scope of the present invention 

[0162] The present invention is not limited in scope to single^ensor or dual-sensor pacemakers, and is not limited 
to IMD's comprlsjig activity orpressure sensors only. Nor is the present inventtam limited in scope to single-chamber 
pacemakers, single-chamber leads for pacemakers or single-sensor or dual-sensor leads for pacemakers Thus, var- 
ious embodiments of the present invention may be practiced in conjunction with more than two leads or with multiple- 
chamber pacemakerSj for example At least some embodiments of the present invention may be applied equally well 
inthe contexts of single^, dual-, triple- or quadruple- chamber pacemakers or other types of IMD's. See, for example, 
U.S. Patent No 5,800,465 to Thompson et al., hereby incprpprated by reference herein in its entirety, as are all U S 
Patents referenced therein 

[0163] IMD 10 may also be a pacemaker-cardioverter- defibrillator ("PCD") corresponding toany of numerous com- 
mercially available implantable PCD's Various embodiments of the present invention may be practiced in conjunction 
with PCD's such as those disclosed in U.S Patent No 5,545,186 to Olson et al., U.S Patent No 5,354,316 toKeimel, 
U S. Patent No. 5,31 4,430 to Bardy, U S Patent No. 5, 1 31 ,388 to Pless and U S. Patent No, 4, ,821 ,723 to Baker et 
al., all hereby incorporated by reference herein, each in its respective entirety 

[0164] Figures 15 and 16 illustrate one embodiment of IMD 10 and a corresponding lead set of the present invention, 
where IMD 10 isaPCD. In Figure 15, the ventrtoularleadtakestheform of leads disclosed in U S. Patent Nos 5,099,838 
and 5,31 4,430 to Bardy. and includes an elongated insulativo lead body 1 carrying three concentric coiled conductors 
separated Irom. ode another by tufatjiar insulative sheaths Located adjacent the distal end pj lead I are ring electrode 
2, extendable helix electrode 3 mounted retractably within Insulative electrode head 4 and elongated coil electrode 5 
Each of the electrodes is coupled to one of the colled conductors within lead body 1 Electrodes 2 and 3 are emptoyed 
for cardiac pacing and for sensing ventricular depolarizations. At the proximal end of the lead is bifurcated connector 
6 which carries three electrical connectors, each coupled to one of the coiled conductors, Defibrillation electrode (elon- 
gated coil electrode) 5 may be fabricated from platinum, platinum alloy or other materials known to be useable in im- 
plantable defibrillation electrodes and may be about 5 cm in length, 

[0165] The atrial/SVC lead shown in Figure 15 includes elongated insulative lead body 7 carrying three concentric 
coiled conductors separated from one another by tubular insulative sheaths corresponding to the structure of the ven- 
tricular lead Located adjacent the J-shaped distal end of the lead are ring electrode 9 and extendable helix electrode 
1 3 mounted retraitably within an insulative electrode head 15. Each of the electrodes is coupled to one of the coiled 
conductors within lead body 7 Electrodes 1 3 and 9 are employed for atrial pacing andf or sensing atrial depolarizations 
Etongatecfcoil electrode 1 9 is provided proximal to electrode 9 and coupled to the third conductor within lead body 7 
Electrode 19 preferably is 10 cm in lengfli or greater and is configured to extend from the SVC toward the tricuspid 
valve In one embodiment of the present invention, approximately 5 cm of the right atrium/SVC electrode Is fecated in 
the right atnum with the remaining 5 (sn located in the SVG At the proximal end of the lead is bifurcated connector 17 
carrying three electrical connectors, each coupled to one of the coiled conductors 

[0166] The coronary sinus lead shown in Figure 15 assumes the form of a coronary sinus lead disclosed in the above 
cited '836 patent issued to Bardy, and includes elongated insulative lead body 41 carrying one coiled conductor coupled 
to an elongated coiled defibrillation electrode 22 Electrode 22, illustrated In broken outline in Figure 15, is located 
within the coronary sinus and great vein of the heart, At the proximal end of the lead is connector.plug 23 carrying an 
electncal connector 41 coupled to the coiled conductor The coronary sinus/great vein electrode 22 may be about 5 
cm in length 

[0167] Implantable PCD 10 is shown in Figure 15 in cotnbinatiCMi with leads 1 , 7 and 41 , and lead connector assem- 
blies 23, 17 and 6 inserted into corihector block 12 Optionally, insulation of the outViard facing portion of housing 14 
of PCD 1 0 may be provided using a plastic coating such as patylene orsilfcone rubber, as is employed in some unipolar 
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cardiac pacemakers: The outward facing portiori. however, may be left unins jiated or some other division between 
insulated and uninsulated portions may be employed The uninsulated portion of housing 14 serves as a subcutaneous 
defibrillation electrode to defibrillate eitherthe atria or ventricles. Lead configurations other that those shown in Figure 
15 may be practiced in conjunction with the present invention, such as those shown in U S Patent No. 5.690,686 to 
Min et aL, hereby incorporated. by reference herein in its entirely. . . 

[0168] Figure 16 is a functional schematic diagram of one embodiment of implantable PCD 10 of the present inven- 
tion This diagram should be taken as exemplary of the type of device in which various embodiments of the present 
invention may be embodied, and not as limiting, as it is believed:that the invention may be practiced in a wide vanety 
of device Implementations, including cardioverter and defibrillators which do not provide anti-tachycardia pacing ther- 
apies. ■ 

[0169] IMD 10 IS provided with an electrode system. If the electrode configuration of Figure 15 is employed, the 
correspondence 1o the illustrated electrodes is as fpllo*vs Electrode 25 in Figure 16 includes the uninsulated portion 
of the housing of PCD 10 Electrodes 25, 13, 21 5 are coupled to high voltage output circuit 27, which Includes 
high voitage switches controlled by CV/defib control logic 3D via control bus 31 Switches disposed withiri circuit 27 
determine which electrodes are employed and which electrodes are coupled to the positive and negative terminals of 
the capacrtor bank (which includes capacitors 33 and 35) during delivery of defibriilation pulses 
[0170] Electrodes 2 and 3 are located on or in the ventricle and are coupled to the R-wave amplifier 37, which 
preferably takes the form of an automatic gain ccntrolied amplifier providing an adjustable sensing threshold as a 
function of the measured R-wave amplitude, A signal is generated on R-out line 39 whenever the. signal sensed.between 
electrodes 2 and 3 exceeds the present sensing threshoid 

[0^ 71] Eiectrodes 9 and 1 3 are located on or in theatrium and are coupled tothe P^tfaveafnpl!fier43, which preferably 
also takes the form of an automatic gain controlled amplifier providing an adjustable sensing threshold as a function 
of the measured P-wave amplitude . A signal, is generated on P-out line 45 whenever the signal sensed between elec- 
trodes 9 and 13 exceetfe the present sensing threshold The gerieral operation of R-wave and R-wave amplifiers 37 
and 43 may correspond to that disctosed in U S. Pat. No 5,117,824,. by Keimel et al , Issued Jun: 2, 1992, for "An 
Apparatus for Monitoring Electrical Physk)loglc Signals", hereby incorporated by reference herein in its entirety 
[0172] Switch matrix 47 is used to select which of the available electrodes are coupled to wide band (0 5-200 Hz) 
amplifier 49 for use in digital signal analysis Selection o1 electrodes is controlled by the microprocessor 51 via data/ 
address bus 53, which selections may be varied as desired Signals fronri the electrodes selected for coupling to band- 
pass amplifier (wide band amplifier) 49 are provided to multiplexer 55, and thereafter converted to multi-bit digital 
signals by A/D converter 57, for storage in random access memory 59 under control of direct memory access circuit 
61 ft/licroprdcessor 51 may employ digital signal analysis techniques to characterize the digitized signals stored irt 
random access memory 5« tg recognize and classify the. patient's heart rhythm ertiploying any of the numerous signal 
processing methodologies known to the art For exaniple, TGRT could be. calculated and output signals generated 
accordingly ^ , „ . 

[0173] The remainder of the circuitry is dedicated to the provision of cardiac pacing, cardioversion and def ibnllation 
therapies, and. for purposes of the present invention may correspond to circuitry known to those skilled in the art The 
following exemplary apparatus is disclosed for accomplishing pacing, cardioversion and defibriilation functions Pacer 
timing/control circuitry 63 preferably Includes programmable digital counters which control the basic t\m intervals 
associated with ODD, VVI, DVI, VDD, AAI, DDI and other modes of single and dual chamber pacing well known tothe 
art Circuitry 63 also preferably controls escape inteivals associated with £mti-tachyarrhythmia pacing in boththe atrium 
and the ventricle, employing any anti-tachyarrhythmia pacing therapies known to the art. 

[01 74] Intervals defined by pacing .circuitry 63 include atrial and ventricular pacing escape Intervals, the refractory 
periods during which sensed P-w^ves and R-waves are ineffective to restart timing of the escape intervals and the 
pulse Widths of the pacing pulses 7ns durations of these intervals are determined by microprocessor 51 , in response 
to stored data in memory 59 and are communicated to pacing circuitry 63 via address/data bus 53 Pacer circuitry 63 
also determines the amplitude of the cardiac pacing pulses under control of microprocessor 51 
[0175] During pacing, escape inteival counters within pacer timing/control circuitry 63 are reset upon sensing, of R- 
waves and P-waves as Indicated by a signals on lines 39 and 45, and in accordance with the selected mode of pacing 
on time-out trigger generation of pacing pulses by pacer output circuitry 65a and 65b, which are coupled to electrodes 
9 1 3 2 and 3 Escape interval counters are also reset on generation of pacing pulses and thereby control the basic 
timing of cardiac pacing functions, including antr-tachyarrtiythmia pacing The durations of the intervals defined by 
escape interval timers are determined by microprocessor 51 viadataJladdress bus .53. Thevalue of the count present 
in the escape inten/al counters when reset by sensed R-waves and R^waves may be used to measure the durations 
of R-R inten/als, P-P intervals, P-R inlen/als and R-P intenrals, which meiasurements are stored in memory 59 and 
used to detect the presence of tachyarrhythmias, 

[0176] Microprocessor 51 most preferably dperates as an Interrupt driven device, and is responsive to interrupts 
from pacer timing/control circuitry 63 corresponding tothe occurrence sensed P-waves and R-waves and correspond- 
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Ing to the gsneralion of cardiac pacing pulses, Thos9 interrupts are provided via data/address bus 53. Any necessary 
mafhematica! calculations to be performed by microprocessor 51 and any updating of the vaiues or intervals controlled 
by pacer timmg'control circuitry 63 take place following such interrupts 

[0177] Detection of atrial or ventricular tachyarrhythmias, as employed in the present invention, may correspond to 
6 tachyarrhythmia detection algorithms such as TCFn". Gonventionally tine presence of an atrial or ventricular tachyar- 
rhythmia is confirrried by detecting a sustained series of short R-R or- P-P intervals of an average rate Indicative of 
tachyarrhythmia or an unbroken series of short R-R or P-P inten/ais. The suddenness of onset of the detected high 
rates, the stability of the high rates, and a number of other factors known iri the art may also bs measured at this time. 
Appropriate ventricular tachyan-hythmia detection methodologies measuring such factors are described in U S Pat. 

JO No. 4,726,380 issued^ to Vollmann, U.S. Pat, No. 4,830,005 issued to Pless et al and U S Pat, No 4,830,006 issued 
to Haiuska et al., all incorporated by reference herein, each in its respective entirety An additional set of tachycardia 
recognition methodologiea Is disclosed in the article "Onset and Stability for Vehtricular Tachyarrhythmia Detection in 
aii Implantabiei Pacer-Cardioverter-Defibrillator" by Olson et al. , published in Computers in Cardiotogy. Oct. 7-10, 1 986. 
IEEE Computer Society Press, pages 167-170, also incorporated by reference herein in its entirety. Atrial fibrillation 

IS detection methodologies are disclosed in Published PCT Application Ser No, US92/02829, Publication No 
W092/18198, by Adams et al., and in the article "Automatic Tachycardia Recognition", by Arzbaecher et al„ published 
in PACE, May-June, 1984, pp 541-547, both of Which are Incorporated by reference herein in their entireties, 
[01 78] In the event an atrial or ventricular tachyarrhythmia is detected and an anti-tachyarrhythmia. pacing regimen 
is desired, appropriate timing intervals for controlling generation of anti-tachyarrhythmia pacing therapies are loaded 

20 from microprocessor 51 into the pacer timing and control circuitry 63, to control the operation of the escape interval 
counters therein and to define refractory periods during whichdetection .of R-waves and P-yyaves is ineffective to restart 
the escape iriterval counters . . 

[0179] Alternatively, circuitry for controlling the timmg and generation of anti-tachycardra pacing pulses as described 
in US. Pat. No: 4,577,633, issued to Berkovits et al, on Mar. 25, 1986, U S. Pat! No. 4,880,005, issued to PIqss et al 

2S on Nov, 14. 1989. U.S. Pal No. 4,726,380, issued to Vbllmann.et al. onfeb 23, 1988 and U S, Pat No 4,587:970, 
issued to Holley et al on May 13, 1 986. all of whkih are incorporated herein by reference in their entireties, may also 
be employed 

[0180] In the event that generation of a cardioversion or defibrillation pulse is required, microprocessor 61 may em- 
ploy an escape inten/al counter to control timing of such cardiovnr.'?ion and defibrillation pulses, as well as associated 

30 refractory periods In response to the dctcclion of strial or vontricuiar fihrillafion or tachyarrhythmia requiring 3 cardi- 
oversion pulse, microprocessor 51 activates cardloveisioiVdofibrillation control circuiiry 30, Which Initiatea charging of 
the high voltage capacitors 33 and 35 via charging circuit 69, under Ihe control of high vpltage charging control line 
71 : The voltage on the high voltage capacitors is monitored, via VCAP line 73, which is passed through rhulliplexerSS 
and in response to reachin.g a predetermined value set by micrcprocisssor 51, results in generation of a logic signal 

3S on Cap Full (OF) line 77 to tenriinate charging Thereafter, timing of the delivery of the defibrillation or cardioversion 
pulse is controlled by pacer timing/control circuitry 63 Following delivery of the fibrillation or tachycardia therapy mi- 
crciprocessor 51 returns the device to .q cardiac pacing mode and awaits the next successive interrupt due to pacing 
or the occurrence of a sensed atrial or ventricular depolarization 

[01 81] Several embodiments of appropriate systems for the delivery and synchronization of ventricular cardioversion 

40 and defibrillation pulses and for controlling the timing functions related to them are disclosed in U S Patent No 
5,188,105 to Keimel, U S Pat No 5,269,298 to Adams et al and U S Pat No 4,316,472 to Mlrowski et al , hereby 
incorporated by reference herein, each in its respective entirety. Any known cardioversion or defibrillation pulse control 
circuitry is believed to be usable in conjunction with various embodiments of the present Invention, however For ex- 
ample, circuitry controlling the timing and generation of cardioversion and defibrillatlcri pulses such as that disclosed 

4S In U S Patent Nn 4,384,585 to Zipes, U S. Patent No, 4,949.719 ts Pless et al., oru S Patent No 4,375,817 to Engle 
et al , al! hereby incorporated by reference herein in their entireties, may also bs employed 
[0l82] Continuing to refer to Figure 1 6, delivery, of cardioversion or defibrillation pulses is accomplished by output 
circuit 27 under the control of controlcircurtry 30 via cbntrol bus31 Output circuit 27 determines whether a monophasic 
or biphasic pulse is delivered, the polarity of the electrodes and which electrodes are involved m deliysry of the pulse 

50 Output circuit 27 also includes high voltage switches which control whether electrodes are coupled together during 
delivery of the pulse Alternatively, electrodes intended to be coupled together dunng the pulse may simply be perma- 
nently coupled to one another, either exterior to or irflerior of the device housing, and polarity may similarly be pre-set, 
as in current implanteJale defibrillators. An example of output circuitry for delivery of bipHasb pulse regimens tomultiple 
electrode systems may be found in the above cited patent Issued to Mehra and in U S, Patent No, 4,727,877, hereby 

ss incorporated by reference herein in its entirety 

[01 83] An example of circuitry Which may be used to control delivery of monophasic pulses is disclosed in U S Patent 
No 5.163,427 to Keimel, also incorporated by reference herein in its entirety Outpuf control circuitry similar to that 
disclosed in U S Patent No 4,953,551 to Mehra et al or U S Patent No 4,800,883 to Winstrom, both incorporated 
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by reference herein in their entireties, may also be used in conjunction witli various embodiments of the present in- 
vention to deliver biphasic pulses. 

[0184] Alternatrvely, IMD 1 0 may be an implantable nerve stimulator or muscle stimulator such as that disclosed in 
U S. Patent No. 5,1 99,428 to Obel et ai., U.S f^tent No. 5,207,218 to Carpentier et al or U.S Patent No 5,330,507 
5 toSchwartz, or an implantable monitoring device such as that disclosed in U.S. Patent No. 5,331,966 issued to Bennet 
' et al all of which are hereby incorporated by reference herein, each in its respective entirety The present invention 
is believed to find wide application toany form of implantable electrical device foruse in conjunction with electrical leads 
[0185] As shown in Figures 12 to 16, conventional implantable devices lO have several electrodes that could be 
used for recording depolarisation and repolarisation signals for the purposes of determining TCfTT or other descriptor 

10 value. For example, in the mo w of Figure 13 signals could be recorded between the atrial tip electrode 21 and the 
ventricular tip electrode 29, between the atriai tip electrode 21 and the IMD ("can") TO and between the ventricular tip 
electrode 29 and the IMD 10 Other electrodes could be placed on the IMD 10, or within the heart 8 or remote from 
the heart connected via additional leads. The signals recorded in three or more directions could be processed to cal- 
culate vectors describing the progress of the depolarisation and repolarisaticffi wavefromts Microprocessor 51 could 
he used to transform the vectors into an optimum dimensional space and descriptors such as TCRT and TMD could 
be determined More complicated I MDs, for example, a pacemaker-cardioverter-defibrlllator, would have an additional 
lead or leads from which further data could be extracted. The "InSync" implantable cardb defibrillator (ICD) which is 
sold by Medtronic, could record data between any two points of the IGD, SVC, three electrocfes in the right ventricle, 
the left ventricle electrode and two atrial electrodes It is desirable to extract the wavefrbnt data from existing leads to 

20 avoid further obstructions being introduced into the heart pathways, for example, through the presence of extra leads 
or thicker leads. In addition, the fewer the number of components, the less likely it is that the device will go wrong 
[0186] In another embodiment, additional electrodes could be positioned on the patient corresponding with the or 
some of the standard independent ECG leads Infprrnatipn from these could be fed to the IMD 1 0 for controlling pacing 
in response to autoricjmicchanges or the detection; of ventricular repolarisation. abnormalities In a further embodiment, 

2S the leads of an IMD 10 could be used to provide data which is fed to external monitoring equipment. This might be in 
the form of a possible device which the patient wears so that the patient can be warned when TCRT values are at 
dangerous levels as a resuit of autonomic changes, for example, during exercise 

[0187] In another embodiment, the IMD 10 may record the data with respect to time which is then fed to an external 
data processing system for analysis, for example, during a check up Such data may provide early indications of heart 
30 problems in addition to that being trsaled by the implantable device 10 These might include ischemia and other con- 
diiions v^hcro damage to the heart muscle is caused which affect the direction of the repolarisation wavefrpnt, change 
the autonomic tone or the response of the autononic system. 

[0188] The descriptor TCRT could be calculated and used by itself in such situations or the ECG data could be 
processed furtherto extract additional descriptor values, which together with TCRT may provide a better diagnosis tool 

35 [0189] In orieembodinnent envisaged in the preseht Invention, a descriptor such as TCRT or TMD, or a combination 
of descriptors, may he used to control a drug delivery mechanism This may he by way of an implantable device or 
possibly via external equipment, for example, when the patient is in an intensive care unit. Such drugs may be used 
to reduce clotting of blood (e,g aspnn) or perhaps reduce tho stress levels In the patient. Sympathetic and parasym- 
pathetic agents could also be administered Where no correction of symptoms is observed, alarms could be triggered 

40 and assistance summoned, ; 

[01 90] As mentioned above, TCRT has been found to provide a useful measure of the autonomic tone in th6 patient . 
Nearly all anaesthetic agents create a sympathetic or a parasympathetic effect in the autonpmic systern of a patient 
It wotild therefore be useful for the autonomic response, of the patient to be tested prior to being anethetised TCRT 
could be monitored w*iilQ a patient performs >felsalva manoeuvre (or some other standard test) to check the autonomic 

4S response, prior to being anethetised for surgery 

[01 91] A third case study was carried to investigate whether QT dispersion could represent properly interlead het- 
erogeneity of ventricular repolarisation The results obtained were compared to determine the extent that the non- 
dipolar components of an ECG differ from the dipolar components, i e the residual energy of the T wave was investi- 
gated The study was as follows. 

so [01 92] The concept ofthe so-called QT dispersion has recently attracted significant attention fromtheclinical research 
community Various methods have been proposed to evaluate QT dispersion from the standard 12 lead electrocardi- 
ogram (EGG) by measuring QT inten/als in individual leads Most frequently, the simple range ofthe QT inten/al meas- 
urements is used There are numerous studies indicating the „clinical value of QT dispersion ArnonaQtiers,. Increased 
QT dispersion has been reported to be associated with QT inten/al prolongation due to drugs of known proarrhythmic 

ss properties and to ba less ihcre&sed on drugs with lesser proarrhythmic effects, to predict mortality in general epidemi- 
ological studies to identify patients who are at greater risk after sun/iving acuta myocardial infarction and lo mark 
therapeutic efficacy in the idiopathic long QT syndrome Recently, however, reports have also appeared challenging 
ttie clinical usefulness of QT dispersicai 
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[0193] Since the introduction of ttie concept of QT dispersion, it has been speculated that tlie increased range of QT 
interval measureinents is caused by the regional heterogeneity of the duration of ventricufar repolarisation. It has been 
proposed that different leads of the standard 1 2-lead EGG project the repolarisation signals of different regions of the 
myocardial tissue and that, consequentiy, increased dispersion is a sign of regional differences in the duration of re- 
s polarisation. Indeed, studies comparing the QT dispersion with the dispersion of these duration of monophasic action 
potentials found a ganeral correlation supporting this hypothesisrit was obsen^ed. that QT dispersion is increasingly 
prolonged with increasing differences in the duration of nKXiophasrc action potentials, recorded at different endocardial 
■ sites. 

[0194] At the same time, these studies do not offer direct proof that increased dispersion of the QT intervals in the 

10 standard 1 2-lead ECGs measures directly the same phenomena as the dispersion of durations of monophasic action 
potentials. If an increased heterogeneity exists in the durations of monophasic action potentials, the repolarisation 
sequence is more disturbed, the vectorcardiographic loop of the T wave is more abnormal and the projections of this 
loop Into the standard ECG leads are more complicated than in normal electrocardiograms. 
[01 95] Recertt studies. have shown that a similar value of QT dispersipri is recorded in full 1 2-lead ECGs and in their 

?5 recisnstruction from orthogonal XYZ leads and that QT dispersion values correlate with parameters of the vectorcar- 
diographic T loop morphology It has consequently beien speculated that the different projections of the T wave vector 
onto the different leads of the standard ECG play an essential role, and that the hypothesis of QT dispersion repre- 
senting a. direct measure of the heterogeneity of ventricular repolarisation durations is flawed Such a concept can 
explain even the earlier obsen/ations of the correlation of QT dispersion with the heterogeneity of monophasic action 

20 potehtiaJ durations However, the studies of orthogonally reconstructed 1 2-lead ECGs and of correlations with T loop 
morphotogy only prove that the projections of the T wave vector play an important role in determining QT dispersion 
but they do not prove that a regional heterogeneity of myocardial repolarisation duration is not involved at all , 
[01 96] To solve the problem of whether QT dispersion is, or is not, associated, witti regional heterogeneity of myo- 
cardial repolarisation, potentially in addition to the T wave vector projection, a dii-ecl study was conducted ccsnparing 

25 QT dispersion wAh. electrocardiographic signals that are not attributable to the orthogonal Vector of tiie T wave In 
1 2-lead ECGs obtained from several clinically well defined, populations, QT dispersion and, using special signal 
processing techniques, the extent of non-dipolar components: by which the individuals leads of a 1 2-lead ECG differ 
from the 3-D vector of the repolarisation signals were measured. 
[0197] The study involved four separate groups of subjects 

30 [0198] The group of normal subjects consisled of 78 normal healthy volunteers (aged 47 + 16 yoars, 23 women) 
wilh normal physical examinalion and normal 1 2-lead ECG At The time of the study, none of the normal subjects was 
on any medication and on the day of the study the subjects were asited to refrain from smelling and from alcohol and 
caffeine intake 

[01 99] The group of hypertrcsphic cardiomyopathy (HCM) patients consisted of 68 patients (mean age 38 ± 1 5 years, 

35 21 women) referred to St George's Hospital London, England, for diagnosis, risk stratification, management of symp- 
toms, and/or follow-up evaluation Following the established guidelines, the diagnosis of HCM was based on the pres- 
ence of left ventricular hypertrophy on 2D echocardiography in the absence of other cardiac or systemic disease that 
may cause left ventricular hypertrophy For ethical reasons, patients were not required to discontinue therapy before 
this study At the time of the' study, '9 patients were on medication with established or potential effects on myocardial 

40 fepolansation (amiodarone, n = 6; and sotalol, n = 5) 

[0200] The group of idlpjaathic dilated cardiomyopathy (DCIVl) patients consisted of 72 patients (mean age 48+ 15 
years, 29 women) The diagnosis of idiopathic DCM was based on enlarged left ventricular diameters (left ventricular 
diameters, diastolic 64 + 10 mm, systolic 51 ± 13 mm) with reduced systolic function without any .underlying causes 
of DCM At the time of the recording, 16 patients were on amlodanane 

4s [0201] Finally, the group of patients with acute myocardial infarction (AMI) consisted of 81 patients (mean age 63± 
12 years, 20 women) Diagnosis of acute myocardial Infarction vuas based on previously published criteria, i e the 
presence of at least 2 of 3 standard signs of (a) typical chest pain, (b) cardiac enzymes elevated more than iwice above 
to normal levels of our laboratoiy, and {c) typjcal ECG changes. History of previous myocardial infarction was recorded 
in 1 3 patients, 44 patients had an anterior irrfarction, and at the time of hospital admission, 68 patients received throm- 

so bolytic therapy At the time of the study,, none of the patients was on an antiarrhythmic therapy 80 patients were receiving 
aspirin, 43 beta-blocl<ers, 32 diuretics, and 34 ACE inhibitors. 

[0202] In the patient groups, subjects were not eligible for this study if in atrial tibriliafion or other non-sinus rhythm, 
in the presence of atrioventricular conduction block, or-with a QRS duration > 120 ms 

[0203] In each subject of each group, 10 serial 1 2-lead ECGs were recorded in a supine resting position using a 
SB digital 1 2-lead electrocardiograph MAC VU by Marquette Medical Systems (Milwaukee, Wisconsin) Each ECG re- 
corded simultaneously all 12-leads for 10 seconds and the serial ECGs were performed one after another without 
removing the electrodes In all subjects, all 10 ECGs were recorded within less than 3 min 

[0204] All ECGs were recorded after careful skin preparation The healthy subjects were recorded after being in- 
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structed to refrain from smoking and caffeine intake on the day of the study. HCM and DCM patients were recorded at 
the lime of presentation at a specialised ouf-patient clinic of our Hospital. The AlVII patients were recorded on day 1 
following the Index infarction 

[0205] Each ECG was stored on a floppy disc (500 Hz sampling at 1 2 bit resolution) and transferred to a dedicated 

s workstation equipped with the QT Guard package (Marquette Medical Systems) which was used to construct the so- 
called median beat of each lead of each eiectrocardiogram. These median beats represent an ideal QRST complex 
of each lead of the ECG and, compared to the native ECG signal, have an improved signal to noise ratio. The median 
beats were further used to measure the QT dispersion and the non-dijDolar components of each ECG. 
[0206] In each eleGtrocardiogram, QT dispersion was measured using the QT Guard package. In each lead of each 

to ECG. the noise of the isoelectric line was measured and compared with the voltage of the peak of the T wave. If the 
standard deviation of the T-P segment signal did not exceed 70% of the maximum t wave amplitude and if the T wave 
amplitude was > 60mV. the lead was measured, otherwise it was excluded from the measurement In each lead, the 
end of T wave was firstly determined automatically using the intersection of the isoelectric line with the tangent to the 
inflection point of the.descending part of the T wave. (The tangent was calculated using least square fit to the 3 samples 

15 above and 3 samples below the inflection point ) These autcsTfjatic measu rements were subjected to visual checks by 
an experienced and electrocardiographlcally trained operator who corrected the automatic measurements manually 
when necessary . During the editing process, the operator was l^ept blindGd with respect to the association of the ECGs 
,«-with the individual groups in the study though clinical ECG ciagnosis was possible in a number of cases. 
;.i[0207] An ECG was accepted for the QT dispersion measurement if the T wave offset was measured in at least 9 

20 .of the 1 2 standard leads. In such a case, QT dispersion was expressed using three different methods, as the range of 
. the QT interval durations in ali measured leads (that is the difference between the maximum and minimum QT Interval 
measured ■ Qld method 1 ), as the standard deviation of the QT interval duratKins in all measurable leads (QTd method 
2), and as the difference between the upper and lower quartjie of the QT iriteryal durations in all measurable leads 
(QTd method 3) Methods 2 and 3 for expressing QT dispersion were used in an attempt tp overcome the technical 

2s problems asaiciated with the simplemeasurementof QT interval range In addition tothese pleasures of Qt dispersion, 
maximum QT interval was taken as the maximum of the QT intervals of all measurable leads Heart rate was also 
derived from each ECG 

[0208] For each method of QT dispersion measurement, Ihs results obtained in the serial electrocardiograms of the 
same subject were averaged and the mean value was used as the true measure ot QT dispersion for the given individual 

30 The representative values of maximum QT interval and of heart rate were obtained for each subject in the same way 
[0209] The concept attributing QT dispersion to the regional differences of myocardial i-epolarisation assumes that 
in addition to the glofaalT wave vector, each lead of the 1 2-lead ECG records signals from a region of the heart (nearest 
to the electrode in case of the precordial leads) which. are. not recorded by any other lead. To quantify the presence of 
such signals, we have measured the non-dipolar compoiients of the 1 2-lead ECG, r e the extent of the residuum of 

35 the T wave which exists after subtracting the T wave veetor 

[0210] More specifically, using the technique described above, in relation to the first study the signals from the eight 
independent leads of the i2-lead ECG (namely leads 1, 11, VI, V2, ,V6) were subjected to the Singular Value Decom- 
position and the electrocarcliogram reconstructed in an orthogonal 8-lead system In such a system, the first lead 
contained the maximum energy in one single direction, the second lead the maximum energy perpendicular to the first 

40 lead, the third lead, the maximum energy perpendicular to the first two leads, etc In this way the energy embedded 
In the first three orthogonal leads corresponded to the energy of the T wave vector while the energy in the remaining 
leads 4 - 8 corresponded to the non-dipolar components summed over all 1 2-leads of the. original ECG (see Figure 2a 
which illustrates the recorded ECG signals and Figure 2b which illustrates those signals wrtienthey are reconstructed 
using Singular Value Decompositlori) TTie Singular V^lue Decomposition is dependent on the region of the ECG signal 

45 for which the optimisatiorr of the orthogonal leads is performed For the purposes of this study, we have optimised the 
Singular \felue Decomposition for the T wave rather than for the QRS complex 

[0211] For each ECG, the propdrtwh between the non-dipolar components In orthogonal leads 4 - 8 and (that is the 
time integral of l^ds 4 - 8 within the T wave) and the energyof the T wave vector in leads 1 - 3 (that is the time integral 
of leads 1 - 3 of the optimised orthogonal system of over the same time) was obtained Similar to the measurement of 
so QT dispersion, these values obtained from the serial ECGs i.i each subject were averaged and the result taken as the 
true measure of the non-dipolar component 

[0212] For the purposes of this study, we term the proportion between the non-dipolar and 3D vector components 
"the relative T wave residuum" 

[021 3] .Subjects were excluded if at least 5 of the 1 0 serial EC(3s eitheir did not provide QT dispersion measu rement 
SB based on the acceptance criteria as abbve, or were rejected by the singular value decomposition package because of 
low signal to noise ratio or other technical reasons 

[021 4] The values of heart rate, Fridericia corrected maximum QT interval (QTc), QT disperskxi and of the relative 
T vrave residuum were compared in individual groups of the study. Since the distribution of the values of the relative 
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T wave residua is not known, the norr-parametric two tail, two sample Marin-Whitney test was used for thiis purpose. 
The correspondence between QT dispersion and T wave residua was examined using Spearman rank correlation 
coefficients which were calculated for the complete st jdy as well as for the individual clinlcaliy defined populations In 
the same way, the con-espondence between T wave residua and heart rate and QTc interval were evaluated 
s [0215] Unless specified otherwise, the data in tables are presented as mean ± standard deviation while in whisker 
charts, data are preeehted as nnean ± standard error of the rheiah, A ja value tessthari 0 05 was considered statistically 
significant. .. 

[021.6] For EGG processing reasons (mainly for T waves of too low amplitude in too many leads), 5 HCM patients 
and 10 DCM patients were excluded from the analysis, in the remaining subjects (78 normal volunteers, 63 HCM 
10 patients, 62 DCM patients, and 81 AMI patients), fewer than all serial ECGs were used in 3 DCM, 1 HCM, and 2 AMI 
patients 

[021 7] Table 7 shows the correlation coefficients between the individual methods for expressing QT dispersion While 
Method 1 (range) was very closely oorretated with Method 2 (standard deviation), Method 3 (inter-quartile difference) 
leads to a somewhat less close correlation although the relationship remains very strong and very statistically significant 
'5 (Figures 17a and 17b) 

[0218] Table 8 shows the differences between the study populations in respect of heart rate and QTc inten/al. Figures 
18a and. 18b shows the values of QT dispersion (Method 1 and Method 3) in the four populations of the study With 
IVIethod 3 (Inter-quartjle range), all the differerices between individual pairs of populations were statistically significant 
or nearly statistically significant. Howevisr, Method 1 (range of QT intervals) and Method 2 (standard deviation) did not 
20 differentiate between normal subjects and DCM patients (p = 0,92 and p = 0 35 for Method 1 and Method 2, respec- 
tively) Note that the results for Method 2 of QT dispersion are not shown in the figure -. although numerically different, 
the values practically reproduced the comparisons with Method 1 
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■ Correlation coefficients between QT dispersion indices 


Group 


QTd 1 vs QTd 2 


QTd 1 ys QTd 3 


QTd 2 VS QTd 3 


NRM ^ 


0 9824 


0 4136 


0 4889 


HCM 


0 9689 


0 6677 


0 7725 


DCM 


09576 


0 6437 


0 7529 


AMI 


09655 


07772 


08685 


Total DODulation 


09812 


0 7212 


0 7942 



QTd. 1 = range of measurable QT intervals,. QJd 2 =: standard deviation of nieasurable QT injen/als, QTd 3 = Inter- 
quartile difterence of measurable QT intervals NRM = norrnal heaflhy volunteers, HCM = hypertrophic; cardiomyopa:fhy 
patients, DCM = Idlopattiic dilated cardiomyopaithy patients, AM| = survivors of acute myocardial infarction 
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Comparisons between heart rate and QTc interval durations in study groups 


Group 


Heart rate (bpm) 


QTc inten/al (ms) 




NRM 


B67±94 


406.5+177 




HCM 


67 6 ±12 9 


447 2 ±26 5 


45 


DCM 


76 0 ±13 0 


429 4 + 35.3 




AMI 


74 1+147 


441 1 + 33 8 



[021 9] All differences between heart rate in individual groups were statistically significant withthe exception of normal 
subjects vs HCM patients, DCM vs AMI patients All differences between QTc iiiten/ai were statistically significant with 
the exception of HCM vs AMI patients NRM = normal healthy volunteers, HCM = hypertrophic cardiomyc^thy patients, 
nCM = idiopathic dilated cardiomyopathy patients, AMI = sunnvors of acute myocardial infarction bpm = beats per 

minute 

[0220] Figure 19 shows the values of the relative T wave residua in the individual populations of the study With the 
exception of the difference between normal subjects and HCM patients (p= 0 14), all the differences between individual 
populations were statistically significant. It should be noted that the values of the relative T wave residuum are very 
small, the mean value in normal healthy subjects being approximately 0 03% which means fftaf in the normal subjects, 
we have found the proportion between non-dipolar and dipolar components of the T wave in the order of 3 in 10,000 
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[0221] The correlations between the relative T wave residua and the measures of QT dispersion is presented in 
Table 9. Corresponding scatter diagrams arc shown in Figures 20a and 20b In individual populations of the study,- 
Statistical significance of the correlation between relative T wave residua and QT dispersion was only reached in HCM 
patients. Note that In the DCM and HCM populations, an opposite relationship between QT dispersbn and relative T 
wave residua was obsen/ed. Statislical significance of the correlation was reached In the total population of the study 
almost certainly because of similarlrertds from norhnal su , ^ 

[0222] Correlation coefficients of the relative T wave residua with heart rate and maximum QTc interval are also 
shown in Table 9. Similar to the rneasuresof QT diversion; the residuurti is related neither to heart rate nortothe QTc 
interval. 

[0223] The findings of this third case study may be summarised as follows: 

a) The non-dipolar components, (i,e electrocardiographic regional heterogeneity) of the repolansation signals are 
rfieasurable in digital 1.2^lead ECGs. 

b) These non-dipoiar componerits. differ in different clinically WBll-defined groups. 

c) The so-caited QT dispersion is unrelated to the non-dipolar components of the T wave Consequently, QT dis- 
persion does n<^ represent a direct measure of regional heterogeneity of ventricular repolarisation 



Table 9 



Correlation coefficients betweein QT dispersion and relativ 


e T wave residuum 


Group 


QTdl vs Twr 


QTd2 vsTwr 


QTdSvsTwr 


HR vs twr 


QTc vs Twr 


NRM 


-0 0446 
NS 


-0 0945 
NS 


-0 0B11 
NS 


0 0794 
NS 


0 2193 
NS 


JICM 


0 2805 
P = 0026 


02882 
P = 0022 


0 3305 
P = 0 008 


-0 2027 
NS 


01322 
NS 


DCM 


-0 1531 
NS 


-0 1755 
NS 


-0 2201 
NS 


■0 0873 
NS 


-0 0317 
NS 


AMI 


0 0771 
NS 


0 0445 
NS 


0 0393 
NS 


01054 
NS 


0.2807 
P = 0011 


Total population 


0 2165 
p = 0 00026 


0 2380 
pi= 6x10-5 


0;2982 
P = 3x10-7 


0.0135 
NS 


0.3270 
P = 2ix10-8 



QTd 1 = range of measurable QT intervals, QTd 2 = standard deviation of measurable QT intervals, QTd 3 - inter- 
quartile difference of measurable QWT intervals HR = heart rate, QTc = Fridericia corrected maximum QT interval 
Twr = relative T wave residuum NRM = normal healthy volunteers, HCM = hypertrophic pardiomyopathy patients, 
DCM = idiopathic dilated cardiomyopathy patients, AMI = survivors of acute myocardial infarction. 
[0224] The findings of this study shows that QT dispersion Is largely caused by the different projebtion of the T wave 
vector in different leads of the standard ECG These findings also show a significant difference between QT dispersion 
in groups with a different morphology of the vectorcardiographic T wave loop Also, this is the first study to show that 
In addition to the projection effects of the T wave vector, no regional components of signal repolarisation play a role in 
determining QT dispersion 

[0225] This case study does not disprove the clinical utility of QT dispersion Indeed, the large number of clmica! 
studies showing the potential of QT dispersion are consistent with the hypothesis thai patients at greater risk (e g 
patients on proarrhythmic therapy or patients with advanced ischaemic heart disease liable to ventricular tachycardia/ 
fibrillation) have a more complex T wave vector and therefore a more complex projection of the T wave vector into the 
individual ECG leads Perhaps, practical aspects of QT dispersion measurement also play a role More complex pat- 
terns of ttie ECG repolarisation signals may clearly lead to increased difficulty with determining the end of the T wave 
and may cmsequently result in the measurement of an increased Ql diversion If this is the case, the more complex 
projections of the T wave vector arc combined with a systematic bias towards increased QT dispersion values in 
patients with disturbed ventricular repolarisation the studies "reporting poor intra- and'inter^ibservef reproducibility 
and poor intrasubject stability of QT dispersion, assessment point in this direction One study found not only increased 
QT dispersion in HCM patients compared to normals but also a lower reproducibilily of Qt dispersion measurement 
in these patients 

[0226] the study, together wnth previously published shows that QT dispersion may relate to nothing more than an 



31 



EP 1 038 498 A2 



expression of T wave loop abnormalltiss, rather than T wave loop morphology. Although technicaify more difficult to 
quantify, the morphology of the T wave loop appears to be a far valuable. ECG factor. 

[0227] I n tfie study it was observed that the relative T wave residuum differs between different clinically well defined 
groups of our study. T wave residyum may be useful clinical^ as a diagnostic tool, either on its own, or more preferably, 
in conjunction with ottier descriptors which could indicate certain tieart conditicsis. It Is possible to take the view that 
tiie relative T wave residuum truly corresponds to the tocaf h 



Claims 

1. A method of characterising ventricular operation of a patient's heart, comprising sensing a plurality of electrical 
signals heart from different spatial positions with respect to the heart during depolarisation and repoiarisation of 
the patient's heart, the plurality of electrical signals monitoring the propagation of depolarisation and repoiarisation 
waves originating in the patient's heart, processingthe plurality of electrical signals to yield a vector which describes 
the propagation direction of one of the depolarisation and repoiarisation waveftorits, and a set of a plurality of 
vectors which describe the propagation direction of the other of the depolarisation and repoiarisation wavefronts 
with respect to time, and determining the vector deviation between the depolarisation and repoiarisation wavefronts 
by measuring the angle between pairs of respective vectors for all combinations of depolarisation vector to repo- 
iarisation vector between predetermined time limits. 

2. A method as claimed in claim 1 , vrfierein the cosine of the angle between- each depolarisation / repoiarisation 
vector pair is calculated. 

3. A method as claimed in claim 1 or 2, wherein anglesare measured between a vector for the repoiarisation wave, 
which cdrfe^ponds to a direction of rnaximum energy of the repoiarisation wave, and each vector from a set of 
vectors describing, the depdarisatioh wave at different time .instances, the set of vectors corresponding to sub- 
stantially the vyhple duration of depolarisation of the patient's heart 

4. A method as claimed in claim 1, 2 or 3, wherein the mean cosine of the angle is calculated for all depolarisation / 
rfipolarisation vector pairs 

5. A method of characterising ventricular operation of a patient's heart, comprising sensing the propagation of depo- 
lansation.and repoiarisation Waves originating in the heatt, determining vectors which arc representative of the 
direction of tiie wavefronts of the depolarisation and repoiarisation waves, and determining the vector deviation 
between the depolarisation and repoiarisation vector? by determining the cosine of the angle between the vectors 
describing the depolarisation and repolarisaition wavefronts 

A method as claimed in claim 5, wherein the vector deviation is a function of. 

a) the cosine of the angle belv.'een hvo vectors, each vector describing one of the depolarisation and repoiari- 
sation wavefronts, 

b) the cosines of the angles between a vector describing either the depolarisation or repoiarisation wavefront 
and a set of vectors describing the other of the depolarisation or repoiarisation wa:vef ront for a plurality of time 
instances, or 

c) the cosines of the angles between a set of vectors describing the depolarisation wavefront for a plurality of 
time instances and a set of vectors descflbing the repoiarisation wavefront for a plurality of time Instances 

A method as claimed in claim. 5 or 6, wherein data is produced for the propagation of the depolarisation and 
repoiarisation wave vyith respect to a first set of axes, tlie data is transformed to a new set of axes defining an 
optimised orthogonal domain having a first axis aligned with a direction of maximum energy for depolarisation or 
repolansation and wherein the vector deviation is measured in the qatimised orthogonal domain 

A method of characterising ventricular operation of a patient's heart, comprising sensing a plurality of electrical 
signals from different spatial positions with respect to the heart during dapolarisafion and repoiarisation of the 
heart, the plurality of electrical signals being sensed by an implantable medicai device and being associated with 
the propagation of depolarisation and repoiarisation waves originating frwrr a patient's heart, processing the elec- 
trical signals to yield a plurality of vectors wrfijch describe the propagatbn direction of a wavefront for a depolari- 
sation wave and a plurality of vectors which describe the propagation direction of a wavefront for a repolarisatbn 
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wave, wherein ventricular operation is characterisad in terms of the cosine of the angle between the plurality of 
vectors for the depolarisation and repoiarisation waves. 

9 A method as clajmed in claim 6, v*erejn at least one of the plurality of vectors describes the propagation of the 
wavefront as a function of time and the mean of the cosine of the angle between pairs of vectors is determined. 

10. An implantable medical device comprising a plurality of medical electrical leads, the leads havrig electrodes for 
sensing electrical signals from different spatial positions iri, on. or near a patient's heart, wherein the device proc- 
esses the electrical signals to yield directions of propagation for depolarisation ahd repolafisatipn waves of a pa- 
tient's heart, calculates the angle of deviation between the depolarisation and repoiarisation waves, and generates 
an output signal corresponding to the angle of deviation 

11. A device as claimed iri claim 1 0, wherein the output signal vanes in accordance with the cosine of the angle of 
deviation 

12. A device as claimed in claim 11, wherein the angle of deviation between the depolansation and repoiarisation 
. . " waves is calculated with respect to time and the output signal varies in accordance with the mean of the cosine 

i of the angle of deviation for a plurality of time instances. 

13. A device as claimed in claim 12, wherein the output Signal varies in accordance with the mean of a set of cosine 
. . vafues for the angle of deviation between the:direcJion of the repoiarisation wavefront for a maximum energy value 

and the di rection of the depolarisation wavefrc«it as a Junction of time for all time instances between start and finish 
points of depolarisation 

14. A device as claimed In any of claims 10 to 13, wherein the output signal is used for at least.pne.of. 

a) control pacing of the patient's heart, 

b) monitor the condition of the patient's heart, 

c) monitor the progression of disease in the patient's heart, 

d) raise an alarm when the angle is outside predetermined limits, 

e) control a drug dispensing pump, and /or 

f) rrwnitor the response of the patient's autonomic system. 

15. A device as claimed in any of claims 11 to 14. wherein the device is selected from the group of. pacemaker, cardi- 
overter, defibrillator, pacemaker-cardioverter-dafibrillator, heart monitor and drug dispensing pump 

16. A device as claimed in any of claims 11 to 15, wherein additional electrical signals are sensed by; 

a) at least one electrode provided on a housing of the device^ 

b) at least one subcutaneous electrdde, and / or 

c) at least one external electrode that is applied to a . patient's body 

17. A device as claimed in claim 16, wrfierein the at least one electrode is coupled to the device by a medical electrical 
lead, by an electncal connection or by a radio frequency transmitter 

18. A method of determining whether a patient is suffering from heart failure, and preferably whether the patient is 
sufiering from erne of hypertnDpic cardiomyopathy; idiopathic dilated cardiomyopathy and acute myocardial infarc- 
tion, comprising measuring the cosine of the angle between propagation directions of depolarisation and repoiari- 
sation waves produced by the patient's heart 

19. A method of determining whether a patient is suffering from heart failure, and preferably whether the patient is 
suffering from one of hypertropic cardiomyopathy, idiopathic dilated cardiomyopathy and acute myocardial infarc- 
tion, comprising sensing a plurality of electrical signals from, different spalialpositions in, on or oaar a patient's 
heart associated with depolarisation and repoiarisation of the patient's heart by employing an implantable medical 
device which is capable of sensing depolarisation and repoiarisation waves and processing the plurality of signals 
to determine the angle between propagation directions of depolarisation and repoiarisation waves produced by 
the patient's heart 
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20. A melhod of monitoring changes in the autonomic system of a patient comprising sensing a plurality of electrical 
signals associated with the propagalron direction of depcdarisation and repolarisatlon waves oi-iginating in the 
patient's heart, determining the cosine of the angle of deviation between the depolarlsation and repolarisation 
waves and monitoring changes in the cosine of the angle of deviation, wherein preferably Scud metfiod is accom- 

5 - plished by employing an implantable medical device virtiich is capable of sensing the waves and processing same. 

21. A method of cJiaracterising ventricular operation, comprising sensing a plurality of electrical signals associated 
with the propagation of a repolarisation. wave originating in the patient's heart, the plurality of electrical signals 
being sensed from different spatial positions on, in or near the patient's heart, processing the plurality of electrical 

10: signals toyleld a plurality of vectors that ai:e representative of the wavefront of the repolarisation wave, and de- 
termining a measure of the spatial variation of the repolarisation wavefront 

22. A method as claimed in claim 21 , wherein the spatial variation is calculated by determining vector contributions 
for the repolarisation wavefront in each of a set of predetermined directions and measuring the angle between 

'5 pairs of vector contributiOTis 

23. A method as. claimed in claim 22, wherein the pre-determined directions, correspond with at least three of the 
standard ECG channels of \, II, VI, V2, V3, V4, VS and V6, and preferably the vector contribution for the ECG 

channel of VI is ignored in the calculation 

20 

24. A method of characterising ventricular operation of a patient's heart, comprising sensing a plurality of electrical 
signals to monitor repolarisation of the heart from different spatial positions with respect to the patient's heart, 
processing the plurality of signals to yield a vector describing the propagation of a repolarisation wave through the 
heart, projecting the vector onto a set of axes to determine vector contributions of the signal vector in the directions 
of the axes, and nieasuring the angle between pairs of vector contributions. 

25. A method as claimed in claim 24, wherein the vector corresponds to a direction of maximum energy of the repo- 
lansation wave 

30 26 A method of characterising ventricular operation of a patient's heart comprising sensing a plurality of electrical 
s'ignals to monitor propagation of repolarisation through the heart from different spatial positions Vvfilh respect to 
the patient's heart, processing the plurality of electrical signals to yield a vector describing the propagation of a 
repolarisation wave with respect to time and with respect to a first set of axes defining an optimum domain space, 
mapping the path of a tip of the vector in the optimum domain space to generate a T-wave loop and calculating a 

35 parameter descnbing the nrrarphology variation of that loop. 

27. A method as claimed In cjalm 26, wherein thd parameter is determined by projecting the T-wave loop on to recon- 
struction vectors corresponding to electrode positions to generate vector contributions in those electrode directiOTs, 
and determining the angle between all pairs of vector contributions 

28. A riiethod as claimed in claim 27, wherein the reconstruction vectors correspond to the position of at least three 
of the standard ECG leads of I, li, VI, V2, V3, V4, V5 and V6, and preferably the vector contribution from the 
standard ECG lead of VI is ignored in the calculation 

45 29. A method as claimed in any of claims 26 to 28, wherein the T-wave loop is rnapped in two orthogonal dimensions 
and the energy of the T-wave loop In thetwo orthogonal dimensions is equalised prior to calculating said parameter 

30. A method of characterising ventricular operation of a patient's heart, comprising sensing a plurality of electrical 
signals associated with the propagation of a repolarisation wave originating in the patient's heart from different 

so spatial positions with respect to the patient's heart, processing data obtained from the plurality of electrical signals 

using a dacompositon technique to yield a plurality of vectors describing the repolarisation wave which are defined 
with respect to a set of axes in an optimised orthogonal domain space, one of the axes being aligned with respect 
to a direction of maximum energy for the repolarisation wave, wherein ventncular' operation is characterised by 
calculating the proportion of repolansation energy not contained within a three dimensional space represented by 

55 the three most significant energy directions 

31. A method of characterising ventncular operation of a patient's heart comprising sensing a plurality of electncal 
signals associated with the propagation of a repolarisation wave originating in the patient's heart from different 
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spatial positions with respect to tha patient's heart, processing tiie plurality of electrical signals to yield a vector 
which is represantative of the wavefront of the repolarisation wave with respect to a first set of axes, transforming 
the vector to a second set of axes defining an optimised orthogonal domain having a first axis aligned with a 
direction ot maximum energy, the domain comprising three dimensrans representing the dipolar cranponents of 
5 the repolarisation wavef rmt vector and at least one further dimension representing the non^dipoiar components 

of the repolarisation wavefronf vector and determining the energy of the non-dipolar components 

32. A method as claimed in claim 31 , wherein said optimised orthogonal domain has eight dimensions and the trans- 
formed repolarisation vravefront vector S has eight corrsponentss,, to SgCprrespo^^^^ dimension, 

to : wherein the vector components are ranked in order.of most.stgnificance with respeict to energy and the non-dipplar 
components are represented by the fourth to eighth componenfe 84 to sg 

33. A method as claimed in claims 31 or 32, wherein the energy of the non-dipcter components es determined for a 
portion of the repolarisation wave corresponding tij a particular region of the heart muscle 

15 

34. A method of determining whether a patient is suffering from heart failure, and preferably a method of w/hether the 
patient Is suffering frotn. one of hypertropic cardiomyopathy, idiopathic dilated cardiomyopathy and acute myocar- 

: dial Inf arctipn , comprising measu ring the energy of the non-dipolar components of a repolarisation wave produced 
by a patient's heart, wherein preferably this is accomplished by employing an implantable medical device which 
so is capable of sensing a repolarisation wave and processing same. 

35. A method of determining depolarisation start and end points for measuring characteristics of a signal representing 
changes in energy during depolarisation of a patient's heart, comprising finding a first peak in the energy of the 
signal corresponding to depolarisation of the patient's heart, determining a point in time.lRp.. corresponding to the 

2S peak energy and determining the maximum energy E^mBx °^ signal at that point, determining a point in time 

t'f,g before tpp and a point in time fpg after ti=!p where the energy of the signal drops to a predetermined percentage 
of the maximum energy, detehnining the depolarisation start point by subtracting a first predetermined, time ihtenral 
from time fpg and detemnining the depolarisation end point by adding a second predatemriined time Interval to 
timet'f^g 

30 

36. A method as claimed in claim 35, wherein sad predetermined percentage of the maximum energy Ep,„ax m 
range ot 50 to 90% of Ef^max' preferably in the range of 50 to 80% of Er^,;, and most preferably 70% of ER^ax 

37. A method as claimed in claim 36, wherein said first and second pfedefermined time jntervals.are in the range of 
3S 38 to 58 msec, preferably in the range of 43lo 53 msec, and most preferably 48 msec. 

38. A method of delemiining repolarisation start and end points for a signal representing, changes in energy during 
depolarisation and repolarisation of a patient's heart, Qomprising: 

40 finding afrfstpeakin the energy of the signal corresponding to depolarisation of the patient's heart, determining 

the maximum energy Er^^^ of the signal at the peak and determining a point in time t-R^ where the energy of 
signal has dropped to a predetennined percentage of Bp^^^^, finding the next peak in the signal energy cor- 
responding to repolarisation and determining the point in time typ where tfiat peak occurs, determining the 
repolarisation start point as a predetermined fraction of the time Interval between 1're and t-rp, 

45 determining the repolarisation end point by determining a vector S2d(U which describes the repolarisation 

wavefront as projected on to a plane spanned by two orthogcmal vectors Ut and Ug: which represent the max- 
imum energy and next nnost energy of the repcdarlsation wave in two orthogonal directions for t, > tjs, the 
vector having a tip which defines a path on said plane, dividing the area defined by the path of the tip of s^o 
(ti) in the plane of and Ug Into a plurality of equa! rectangular cells, assigning a measure D, to each cell 

so dependent on the time spent by thetipof S2D(ti) in the i"^ cell, discarding cells having the measures Dj = 0 and 

ordering ail other cells in respect of D,, determining a threshold value 0,^ of D, which is greater than the mean 
value of D,, and determining the end point of repolarisation tjE as a point at which D| > D,f, 

39. A method as claimed in claim 38, wherein sad predetermined fraction of the time intenifal between t'f,E and t^p is 
ss in the range ]4 to Yz preferably said predetermined fraction of the time interval between and \jp is 

40. A method as claimed in any of claims 1. to 9, or 20 to 33, wherein the plurality of electl:ical signals are measured 
from electrodes positioned in different spatial positions corresponding to the standard ECG leads of i, II, VI, V2, 
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V3. V4, V5 and VS. 

41 . A method as claimed in any of claims 1 to 9, 20 to 33 or 40, wrfierein data from the electrical signals is decomposed 
into an optimum domain space using singular value decomposition 



42. A method of detecting whether a patient is healthy or sick, wherein more than one parameter characterising ven- 
tricular operation as claimed. in any of clainis 1 to 9, 21 to 33, 40 or 41 is calculated 

43. A method of detecting abnonrnalities of ventricular fepblarisation in a patient using a method of characterising 
ventricular operation asclaimed in ariy of claims 1 to 9, 21 to 33, .40or41,.or.a method of monitonng changes in 
the autonomic system of the patient as claimed in claim 20, .wherein preferably said method is a method of detecting 

the onset of ischemia 

44. A method of determining whether a patient is suffering from heart failure, and preferably a method of determining 
whether the patient is suffering from one of hypertropic cardiomyc^athy, idiopathic dilated cardiomyopathy and 
acute myocardial Infarction, comprising measuring a parameter characterising ventricular operation as claimed in 
any of claims l to 9, 21 to 33, 40 or 41, wherein preferably the method is accOTnplished by employing an implantable 
medical device Which is capable of sensing depolarisation iand.repolarisation waves, and processing same 

45. A method of categorising subjects using a method of characterising ventricular operation as claimed in any of 
claims 1 to 9, 21 to 33, 40 or 41, or using a method of monitoring changes in the autonomic system of a patient 

as claimed in . claim 20 

46. A method of mon itoring the autonomic system of a subject, whereiria method of characterising ventricular operation 
as claimed in any pf claims 1 to 9, 21 to 33, 40 or 41 , is ussd to measure changes in the autonomic system of the 
subject, wherein preferably, the method is accomplished by employing an implantable, medical device which |e 
capable of sensing depolarisation and rspolarisation waves, anci processing same. 

47 A method as claimed in claim 46, wherein changes in the autonomic system of the patient are measured to. 

a) monitor Ihe progress of a disease in a subject, 

b) monitor the influence of drugs on the autonomic system of a subject, 

c) controt the rate of pacing for a pacemaker, 

d) control a drug dispensing pump and / or 

e) test the response of the autcsiomic system of a subject. 

4B. A method as claimed in claim 47, wherein changes in the autonomic system of the patient are measured while the 
patient executes a pre-determined procedure to effect a change in autonomic system, wherein preferably the 
autonomic system is measured while Valsalva manoeuvre isperformed by the patisnt.and/or while aset of postural 

changes are executed by the patient 

49. A method of detecting whether a patient is healthy or sick comprising sensing a plurality of electfica;! signals from 
different spatial positions in, on or near the patient's heart during depolarisation and repolarisation of the patient's 
heart and ctetemiining a parameter. 

a) Which measures vector deviation of the depolarisation and. repolarisation waves in terms of Total Cosine 
R_to_T (TCFIT), 

b) whkii measures the spatial variation of a repolarisation wavefronl in terms of T-wave Morphology Dispersion 
(TMD), and/or 

c) v**ilch measures local repolarisation defects as a function of the energy of the non^dipolar energy compo- 
nents of a repolarisation wave. 

50. A system for measunhg characteristics of ventricular qDeratfon of a patient's heart in accordance with any of the 
methods of claims 1 to 9, 21 to 33, 40 or 41 , or a system for monitoring changes in the autonomic system of a 
patient as claimed in claim 20, comprising a plurality of electrodes for sensing a plurality of electrical signals in 
different sjjatial positions in, on or near to a patient's heart to monitor the propagation of depolarisation and repo- 
larisation waves originating In the patient's heart, a microprocessor for processing the plurality of electrical signals 
to detemiine a characteristic of ventricular operation as claimed in any of claims 1 to 9. 21 to 33, 40 or 41 , or 
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condition of the patient's autonomic systGm as claimed in claim 20, and an indicator to convey the value of tlie 
mGasured characteristic, or condition of tie patient's autonomic system, to an observer, wherein preferably the 
Indicator is a visual display and/or an audible alaim 

51. A system as claimed in claim 50, wherein the system senses the plurality of electrical signals via. 

a) at least one electrode attached to a medical electrical lead of an implantable medical device selected from 
the group of pacemaker, cardioverter, defibrillator, pacKnaker^rdioverter-defibrillator and heart monitor 
which IS implanted within the subject's body; 

b) at least one electrode provided on a housing of an implantable medical device, 

c) at least one su'bcufaheou^^'electrbd~e:"ah^ " ■ ' 

d) at least one external electrode applied to the patient's body 

52. A system as claimed in any of claims 50 or 51 , wherein the system has eight electrodes which correspond to the 
standard electrode positions of I, II, VI, V2, V3, V4, V5 and V6 

53. An implantable medical device which measures charactenstics of ventricular opeiration of a patient's heart in ac- 
cordance with any of the methods of claims 1 to 9, 21 to 33, 40 or 41 , or monitors changes in the patient's autonomic 
system as claimed in claim 20, comprising a plurality of electrodes for sensing a plurality of electrical signals in 
different spatial positions in, on or near to a patient's heart to monitor the propagation of depolarisation and repo- 
larisation waves originating in the patient's heart, and a microprocessor for processing the plurality of electrical 
signals to measure a characteristic of vantricular operation as dainried in any of claims 1 to 9, 21 to 33, 40 or 41 , 
or "condition of the patient's autonomic jsysterri as claimed in claim 20, and to generate an output signal corre- 
sponding to. said rtieasured characteristic 

B4. Apparatus as claimed In claim 53, wherein the apparatus Is at le&st one of: 

a) a pacemaker or pacemaker-cardioverter-defifarillator having a pacing rate which is controlled in response 
to the calculated characteristic of ventricular operation,' 

b) a cardioverter or a defibrillator, the Operation of which is controlled In response to the calculated characler- 
istic ol ventncular operation, 

c) a monitor for detecting abnormalities in ventricular repoiansation, the device preferably being able to trigger 
a warning signal when the calculated characteristic of ventricular operation is within certain limits, 

d) a monitor for measuring changes in the autonomic system of a sjbject, the device preferably being able to 
trigger a warning signal when the calculated characteristic of ventricular operatibn is within certain limits, and 
/or 

e) a device for controlling a drug dispensing pump 

55. A computer program, compnsing software code portions for performing the method as claimesd m any of claims 1 
to 9, 20 to 33, 40 or 41, a computer program product which is directly loadable into the intenial memory of a digital 
computer, comprising software code portions tor performing the method as claimed in any of claims 1 to S, 20 to 
33, 40 or 41 , a microprocessor which is programmed with software code portions for performing the method as 
claimed in any of claims 1 to 9, 20 to 33, 40 or 41 , or a. digital computer which Is progrannmed wfth software code 
portions for performing the method as claimed In any of claims 1 to 9, 20 to 33, 40 or 41 
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The Search Division considers that the present BuropBan patent appiicatiori doss not comply with the 
requirements of unity of invention and relates to several Inventions or groups of inventions, namely: 

1. Claims: 1-17,20 atid 48,41.45,46,48.55 (part) 

A method of and device for characterising ventricular 
operation comprising; measuring EC6 signals, yield vectors 
that describe propagation of depolarisation and 
repoTarisation wavefront with respect to time and measuring 
angles between pairs of respective vectors. 



2. Claims: 21-25 and 4G, 41. 45,46,48,55 (part) 

A method of characterising ventricular operation comprising: 
measuring ECG signals, yield vectors that describe 
propagation of repolarjsation wavefront and detemining a 
measure of the spatial variation of the repolarisation 
wavefront;/. 



3. Claims: 26-33 and 40,41,45,46,48,55 (part) 

A method of characterising ventricular operation comprising: 
measuring EC6 signals, yield a vector that describe 
propagation of repolairisation wavefront with respect to 
time, generate a T wave Loop and calculate parameters 
therefrom. 



4. Claims: 35-37 and 40,41,45,45,48,55 [part) 

A method of determining depolarisation start and end points 
COTiprising: finding a peak in the energy signal 
corresponding to maxiinupi energy and defining the start and 
end points as predeteimi ned percentages thereof. 



5. Claims: 38, 39 and 40.41,45,46,48,55 (part) 

A nstbod of deteirmining repolarisation start and end points 
comprising: defining the start point as a fraction of the 
time difference between two peaks in the depolarisation 
energy and defining the end point from a vector describing 
the repoTarisatTOTi wave in. two orthogonal directions. 



6. Claims: 59-54 and 40,41,45,46,48,55 (part) 

A system and device for measuring characteristics of 
ventricular operation comprising: a plurality of electrodes, 
a microprocessor, an indicator to convey measured values, 
said indicator being a display and/or an alarm 
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